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ABSTRACT 

This thesis uses Flightlab to model and analyze the flight dynamics of the SH-60B 
Seahawk helicopter. Flightlab runs on computers utilizing the UNIX™ operating system and 
is used for design, analysis and testing of an aircraft using non-linear modeling techniques. 
It is capable of modeling conventional main rotor-tail rotor and tandem rotor helicopters and 
tilt rotor aircraft. A procedural guide for modeling and analyzing a single main rotor 
helicopter is presented using the SH-60B. The non-linear response from lateral and 
longitudinal cyclic, main rotor collective and tail rotor collective inputs are presented. 
Flightlab is also capable of reducing the non-linear model to a linear model. The linear and 
non-linear models are then compared. The purpose of this thesis is to present a guide for 
using Flightlab to model and analyze an existing helicopter design, and also to have in place 
a well tested model to be used for further research. 
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I. INTRODUCTION 


A.BACKGROUND 

The advent of high speed processors and realistic computer generated dynamic 
displays have made the computer an indispensable tool for rotorcraft design and analysis. 
The hardware though is only one part of the equation, the other being software. The four 
major helicopter manufacturers in the United States, McDonnell Douglas, Sikorsky Aircraft, 
Bell Helicopter Textron and Boeing Helicopters have all invested heavily in computer flight 
simulation, modeling, and engineering analysis. The most ambitious is probably Sikorsky 
Aircraft's expansion of its simulation facility to 300,000 square feet, mainly to support the 
LHX program [Ref. 1: p. 153], The selection of the Boeing/Sikorsky team to build the new 
Light Helicopter (LHX) for The U. S. Army was the result of competition based solely on 
the flight simulation of an engineering design [Ref. 2: p. 1]. This was accomplished prior 
to the first prototype being built. High end rotorcraft modeling software is used across all 
lines of analysis and design, not just flight modeling. It is also used for hardware in the loop 
testing, piloted simulations and rapid prototyping of new system designs. 

The large interest in computer flight analysis and simulation has led to the 
development of a software package called Flightlab. This program provides greater fidelity 
in modeling rotorcraft over purely mathematical based models usually used in MATLAB™ 
or other computational software. Flightlab’s strength in this area is due to the fact that 
rotorcraft models are constructed using physical components, i.e. a hinge, wing or distributed 
mass to mention a few. The computations to handle the reaction of aerodynamic and 
gravitational forces between the model components is handled by the software. The user 
does not have to derive a mathematical expression to represent the dynamics. This approach 
leads to a nonlinear analysis which is valid over the entire flight envelope as opposed to a 
linear model which is valid only within the small linear region. 

A second benefit is that the rotor system degrees of freedom are calculated in addition 
to the six degrees of freedom of the body. This permits a more detailed analysis to be 
conducted on the model and thus an increase in the accuracy of the results. Flightlab is also 




capable of real time flight simulation for piloted evaluations. This requires the use of parallel 
processors. 

B. OVERVIEW OF FLIGHTLAB 

The Flightlab program, current version 2.5, is written in a higher level language 
which uses a combination of FORTRAN and C computer languages. This fact is transparent 
to the user, but a familiarity with the UNIX™ operating system is required. 

There are three programs within Flightlab for nonlinear modeling of aircraft. Model 
Editor, Gscope, and Model Analysis. The Model Editor is a high level modeling tool that 
allows for quick construction of a model through the use of standardized model structures. 
Each of the eight standardized model structures, also called supercomponents, facilitate for 
a modular design which allow for easy changes in model design. In addition, the Model 
Editor provides for selective fidelity modeling options. For example, the rotor system can 
be modeled as either a rigid or elastic system. Pertinent data for selected supercomponents 
is entered into data templates either directly as numeric data or the names of data files for 
more involved data, for example coefficient of lift versus angle of attack of the fuselage. 

Gscope provides an environment similar to SIMULINK™ for constructing systems 
or components of the aircraft. Once a system is constructed in Gscope, it is linked to the 
model created in the Model Editor. Gscope is most often used when a particular system of 
the aircraft under study differs from the generalized model structure present in the Model 
Editor. It is even possible to construct an entire model in Gscope without using Model 
Editor, although with a much higher workload. 

The Model Analysis program contains numerous scripts accessed in a menu format 
for analyzing the dynamics of the model. It is a very extensive analysis program and can 
also be modified to meet the needs of the user. The second aspect of the Model Analysis is 
a command line interface with many of the basic features of MATLAB™ for matrix 
manipulation and plotting of data. Flightlab uses the interactive plotting capability of 



c. 


MOTIVATION 


The Flightlab program provides for dynamic real time simulation of an aircraft model 
that can be used in a design course within a one semester time frame, but the rapid 
improvements in the program have outpaced the development of user manuals provided by 
the company. The developers of Flightlab provide five user manuals. Scope Users Tutorial 
[Ref. 3]; Component Reference Manual [Ref. 4]; Scope Theory Manual [Ref. 5] ; Scope 
Command Reference Manual [Ref. 6]; and the latest, Flightlab Training Manual [Ref. 7]. 

The first four manuals were written for one of the first versions of Flightlab, which 
only had Gscope as its modeling and analysis tool. It required a great deal of time and 
expertise to implement and analyze a design and thus would not allow for the design to be 
completed in one semester unless the user was already well versed in its use. In addition, a 
previous Naval Postgraduate School thesis. Procedural Guide For Modeling and Analyzing 
the Flight Characteristics of a Helicopter Design Using Flightlab [Ref. 2], covered the use 
of Gscope for the complete modeling of an aircraft, where references 3 thru 6 covered much 
of the basics of modeling various items, but they do not show how to model a complete 
aircraft. 

Flightlab version 2.4 introduced two additional tools. Model Editor and Model 
Analysis. Reference 7 is a tutorial on how to use these two tools. This improvement 
drastically reduced the amount of time needed to implement and begin testing a model. The 
model editor contains generalized aircraft components of varying degrees of fidelity, which 
are chosen based on simulation requirements. Data is then entered in a graphical user 
interface format, the file saved and analysis begun. Flightlab version 2.5, as of now, has no 
manuals describing it's improvements over previous versions. 

The purpose of this thesis is two fold. First, to provide a structured, logically 
sequenced instruction manual for the implementation and testing of an existing aircraft 
design. Also, to explain some of the finer points of the program that are usually experienced 
only after using such a program. 

Second, to model and analyze the flight dynamics of the SH-60B Seahawk helicopter 
and compare the results with those obtained by Sikorsky Aircraft using their analysis 





program, GenHel. Aerodynamic and control system data for the SH-60B was supplied by 
Sikorsky Aircraft. The flight characteristics of the closed and open loop model will also be 
studied. 



II. FLIGHTLAB TUTORIAL 


A. GETTING STARTED WITH FLIGHTLAB 

Flightlab is currently installed on the Naval Postgraduate School Department of 
Aeronautics SGI Indigo II workstation hawkeye, but can be easily accessed by other 
aeronautics workstations utilizing the UNIX™ operating system environment. The 
procedure by which will be explained shortly. Setting up an account to run Flightlab requires 
modifying the .login and .cshrc files as in Table I. 

.login file 

setenv FL_DIR /usr/local/flightlab-2.5 
source SFLJDIR/SOURCE.ME 
setenv PS_PRINTER 'hp3si_l' 
setenv FL_TERM xterm 

.cshrc file 

set path=( /usr/local/flightlab-2.5 /usr/local/ 
flightlab-2.5/bin Spath) 
xhost hawkeye.aa.nps.navy.mil 

Table I. Login Files Modification 

After modifying the .login and .cshrc files, the process of starting Flightlab will vary 
slightly depending on whether the user is physically located at hawkeye or at another 
workstation. If at hawkeye, the program is started by typing "flightscope &" at the users 
home directory. The ampersand is optional, but allows the shell window to be used while 
the program is running in the background. This is convenient as it allows for editing of files 
and viewing the various Flightlab windows at the same time. If at a workstation other than 
hawkeye, follow the procedures in Table II. 





rlogin hawkeye.aa.nps.navy.mil 
enter user name and password 
setenv DISPLAY phantom.aa.nps.navy.milrO.O * 
flightscope & 

♦Enter the terminal name as appropriate, phantom is used as and example. Ensure that 
the :0.0 is included immediately following the terminal name. 

Table II. Remote Login Procedures 

Flightlab will begin by displaying the main window as shown in Figure 1. The three 
main tools of the program, Model Editor, Gscope, and Model Analysis are each accessed by 
clicking once on the appropriate square with the left mouse button. 



Figure 1. Flightlab Main Window 











B. MODEL EDITOR 

The Model Editor is where a majority of the data for the aircraft under study will be 
entered. When selected, it brings up a model tree which contains all model structures 
necessary for the complete modeling of an aircraft. Figure 2 is a compressed version of the 
full model tree, but does show the eight supercomponents and a number of selected 
subcomponents. The eight supercomponents are, EnvironAtmos, Rotorl, Rotor2, Wing, 
Propulsion, Airframe, Flight Control, and Aero Interference. 



Figure 2. Model Editor Tree 
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Construction of a model begins by first selecting the Rotorcraft option at the top of 
the tree with one click of the left mouse button. The Rotorcraft option will then become 
enclosed in a red box indicating that it has been selected and will become part of the model. 
This is true for all supercomponents and subcomponents selected. 

Next, the appropriate supercomponents and subcomponents are selected in the same 
fashion. Flightlab will not allow a subcomponent to be selected prior to it's corresponding 
supercomponent being selected and enclosed in a red box. This is what is refereed to as a 
top down design. Hereafter supercomponents and subcomponents will be generically refered 
to as components. 

After all desired components are selected, a double click with the left mouse button 
on any of the previously selected components will bring up it's corresponding data template, 
although not all components have a data template. A representative data template is shown 
in Figure 3. 



It is important to note that not all of the eight supercomponents must be selected in 
order to start analyzing the model. Flightlab is flexible enough to allow each system to be 
constructed in sequence and then tested to allow for easy debugging both in syntax and logic 




implemented. For example, the EnvironAtmos and Rotor 1 components can be selected, the 
appropriate data entered, and then the model analyzed as an isolated rotor in the free stream. 

Figure 3 shows that the required data is entered either as numerics or the name of a 
file which contains the necessary data in table form. The directory path for data files is 
specified in the “Path” menu on the upper left comer of the model tree. Beneath “Path” 
select the “Path for Tables...” option. Flightlab will respond with a window to enter the 
desired path. All data files refereed to in the model must be present in this path. 

Data table files are called by one of two methods. File names with the .sav suffix are 
read in with the “read” command and stored in RAM to the variable assigned. For example, 
the table corresponding to blade_cg.tab in Figure 3 would be: bchord=read("directory 
path/blade_cg.tab "). Comments in .sav files are delineated by the pound sign, #, directly 
preceding the comment. 

The second method corresponds to files with the .tab suffix. Here the “load” 
command is used and the variable name is included in the data file as the first readable 
element. The syntax for the load command is: load("directory path/cl 1.sav"). Comments 
in .tab files are delineated by the pound sign , #, directly preceding the comment. 

The two methods for calling a data file also differ in the way the data is read from the 
file and entered into the matrix. The load command requires that the variable name for the 
data file be the first readable element in the file, immediately followed by three numbers 
which specify the size matrix the data will be entered into. The first is the number of rows, 
the second the number of columns and the third the depth of the table, usually zero. The data 
file will be read from left to right top to bottom filling in the matrix columns first, i.e. first 
row first column, second row first column to the end of the column and then starts on the 
second column. Table III demonstrates how to enter the data of a two dimensional array into 
a data file to be called by the load command. 

The read command only requires two numbers to specify the matrix size and does not 
require the variable name in the data file, as it is assigned by the calling command. The data 
is read from the file left to right top to bottom as before, but the data is entered into the 
matrix in a different fashion. The rows are filled in first followed by the columns, i.e., first 




row first column, first row second column to the end of the first row and then starts on the 
second row. Table IV demonstrates how to enter the data blade chord vs blade station to be 
called by the read command. 

After all the data is entered, click on the “mark save and done” and then “close” 
boxes on the bottom of the data template. The template will close and its corresponding 
component on the model tree will now be enclosed in a red box with the left half colored 
yellow. This indicates that data entry is complete at this node. If a data file name entered 
on the data template is not found in the directory path specified, Flightlab will give an alert 
stating so, but the alert can be overridden and the node saved regardless. 


Data given in following format 

input 1 = sideslip angle -10 to 10 delta 5 
input 2 = angle of attack -10 to 10 delta 10 
output = pitching moment 

angle of attack = -10 

20.0 18.0 0.0 0.0-7.0 

angle of attack = 0 

-39.0 -5.0 0.0-16.0 -55.0 

angle of attack = +10 
-56.0-15.0 75.0 128.0 171.0 


Data file format 
variable name 5 3 0 

20.0 18.0 0.0 0.0 -7.0 

-39.0 -5.0 0.0 -16.0 -55.0 

-56.0 -15.0 75.0 128.0 171.0 

If the variable name is called up in the GScope command line, the matrix will be 
displayed as follows: 

20.0 -39.0 -56.0 

18.0 -5.0 -15.0 

0.0 0.0 75.0 

0.0 -16.0 128.0 

-7.0 -55.0 171.0 _ 

Table III. Data File Structure for Load Command 
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Data given in following format 

Blade station .041 .046 .100 .200 .300 .400 .500 .600 .700 .800 .900 1.00 


Chord 0.00 1.73 1.73 1.73 1.73 1.73 1.73 1.73 1.73 1.73 1.73 1.73 


Data file format 
12 2 

0.0410 0.00 0.0460 1.73 0.1000 1.73 

0.2000 1.73 0.3000 1.73 0.4000 1.73 

0.5000 1.73 0.6000 1.73 0.7000 1.73 

0.8000 1.73 0.9000 1.73 1.0000 1.73 

If the variable name is called up in the Gscope command line, the matrix will be 
display as follows: 

0.4100 0.00 

0.0460 1.73 

0.1000 1.73 

0.2000 1.73 

0.3000 1.73 

0.4000 1.73 

0.5000 1.73 

0.6000 1.73 

0.7000 1.73 

0.8000 1.73 

0.9000 1.73 

1.0000 1.73 


Table IV. Data File Structure for Read Command 

The model can now be saved by accessing the “File” menu at the upper left corner 
of the model tree and selecting the “Save model as...” option. Prior to saving the model 
though the “Path for Local Database...”, beneath the “Path” option must be set. Flightlab 
comes with a selection of completed ready to run models that are located in a central 
database. The path option for these are: /usr/local/flightlab-2.5/flightscope/models/uh60. 
Generally, models to be saved in the users account will have the following path: /d4/user 
name/flightscope/models/rotorcraft. 
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c. 


GSCOPE 

When selected, the Gscope window will be displayed as shown in Figure 4. The 
purpose of Gscope is to provide a tool to build dynamic models, enter associated data and 
script files, and conduct analysis on the model using the Scope language. Scope consists of 
a set of commands and functions that support linear systems analysis, parameter 
identification, eigenvalue analysis, fourier transforms and general matrix manipulation [Ref. 
1: p. 22], Scope Command Reference, Reference 4, gives the basic operations, syntax, 
functions, and commands supported by Scope. These basic commands can be combined to 
form new functions, script files, that can be edited and then executed. These user defined 
script files can be as simple as a series of repetitive calculations to executing an entire 
rotorcraft model and plotting the desired output. 



Figure 4. Gscope Main Window 


Gscope has five window options, main window, model window, editor window, help 
window, and plot window. All are accessed via the “Windows” pull down menu at the top 
of the main window. 
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1. Main Window 

Referring to Figure 4, the main window has three basic areas, the display, history and 
command entry screens. The display screen is the upper area of the main window below the 
pull down menus. This is where the commands and output of the Gscope program are 
displayed. 

The middle area is the command history screen. It is a buffer which holds the twenty 
most previous commands entered. Previous commands may be copied to the current 
command line by clicking on the desired one with the left mouse button. After clicking on 
a command, the others can be scrolled through by using the up and down arrow keys. The 
command line screen is just below the S> at the bottom of the main window. Single line 
entry of Scope commands are entered here. 

2. Model Window 

The model window can be accessed by selecting it beneath the “Windows” pull down 
menu at the top of the main window. The model window shown in Figure 5 is used to build 
and design the model. 
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Construction of a model is accomplished by selecting and placing components in the 
model screen, connecting them and then entering the pertinent data in dialog boxes. 
Components were designed to model a specific physical entity such as a hinge, a distributed 
mass, or a wing, to mention just a few. For convenience, components are grouped into 
classes. The five classes are kinematic, aero, control, analytic and solution. The icons of the 
components within the five classes are included in Appendix A. Each component class is 
accessed on the left side of the model window. Figure 5 shows the control class. To display 
a different class, place the cursor over the class name and press the left mouse button, a menu 
of the five classes will be displayed. Simply move the cursor to the desired one and release 
the mouse button. 

Components are selected by first placing the cursor over the component in the group 
selection area to the left and clicking the left mouse button. It is then deposited by moving 
the cursor into the model screen and clicking again. Components are connected by placing 
the cursor on the component, holding down the center mouse button and dragging the cursor 
to the next component, at which time the button is released. Entering of data for components 
is accomplished by double clicking on the desired component with the left mouse button. 
A dialog box will be displayed showing the type of component and the mandatory data fields 
to be entered, Figure 6. Generally, each component has more variables associated with it 
than shown in it’s dialog box. The others are computed during execution of the model. 
Reference 2 lists the component variables associated with each component. 



Figure 6. Dialog Box 
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Data assigned to each component variable can be either a numeric or point to a 
variable name by using the “&” operator. Variables assigned to component variables in this 
way are called “scope component variables”. The other type of variable in Flightlab is 
known as a “scope free variable”. The use of scope component variables adds flexibility to 
the model design. Variable values can be changed once in the .prolog file, vice changing it 
in several different components throughout the model. The use of scope component 
variables also negates the need of regenerating a script file each time a minor change is 
made. Refer to Figure 6 for an example of the “&” operator use. 

Editing of components can be done by either selecting the “Edit” pull down menu or 
highlighting the desired component and then pressing the right button which will display a 
pop up menu. A component is highlighted by clicking once on it with the left mouse button. 
If upon selecting this pop up menu the display freezes, ensure that “Num Lock”, “Caps 
Lock”, and “Scroll Lock” options are all de-selected on the keyboard, i.e., the lights above 
these options are off. After turning them off the display will free up. Grouping of 
components is done by selecting the “Group” icon in the upper left comer of the model 
window, moving the cursor into the model screen, hold down the left button and as the cursor 
is dragged a box will be displayed. Adjust the box until it encloses all desired components 
and then release the button. 

The model window has five pull down menus, most of which are self explanatory, 
but several are unique to Flightlab and will be discussed here. 

Beneath the “File”pull down menu is an item called “Generate Script...”. This option 
is chosen after a model is complete and ready to be analyzed. Upon selecting it and entering 
the desired file name, a text file is generated with the suffix .exc. This text file can be 
viewed with the editor. The script file generated cannot be executed by itself to analyze the 
model. It will be called by another script file, generally with a .def suffix generated by the 
model editor, which will generate results. 

The “Group” pull down menu has three options, “Goto Parent Group", “Ungroup 
Selection”, and “Ungroup All”. These options allow the user to display the various levels 
which may exist in the model hierarchy, and also to move selected component(s) or entire 




groups of components to other levels. Model hierarchy will be discussed in detail in Chapter 

in. 

The “Goto Parent Group”, when selected, will bring up the display of the group 
which is one level higher in the model hierarchy. To move in the opposite direction i.e. 
to the display the group one level lower in the model hierarchy, double click on the desired 
group icon. 

The “Ungroup Selection” menu item will simultaneously move the selected, 
highlighted, component to the next higher level and display that level. If it is desired to 
place a component into another existing group, drag it and place it on top of the group icon, 
the component will no longer be displayed when it has been successfully transferred to the 
new group. Double click on the group icon to display it and verify that the component has 
moved. These two features allow a component to be connected to another which is located 
in another group. For example, select a component so it is highlighted and then choose 
“Ungroup Selection” beneath the “Group” pull down menu. Repeat this until the component 
is at the desired level. Now connect the components. Return the component to the original 
group by placing it on top of the group icons, in succession, as described above. 

The “Ungroup All” menu item will ungtoup the present level displayed on the model 
screen and place all components in the next higher level. The group icon will remain, but 
it will be empty, if not needed it can be deleted. 

3. Editor Window 

The editor window provides a basic text editing capability. Scope files can also be 
executed from here if it is currently displayed in the editor window. The execute option is 
below the “File” pull down menu. 

4. Help Window 

Allows for on line help with scope commands, syntax and components. 

5. Plot Window 

Provides a graphical display of data as generated from the main window under scope 
execution. 
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D. MODEL ANALYSIS 

Model Analysis provides an interactive menu driven format for the analysis of a 
model created under Model Editor and Gscope. After a short introductory message, the Main 
Menu will be displayed which is the top level of the analysis hierarchy. The general 
procedure for analyzing a model is a follows: 

• Load the model 

• Set initial conditions 

• Prescribe desired simulation 

• Execute simulation 

• Plot results 

• Save results, if desired 

The Model Analysis program operates similarly to Gscope in that it generates scope 
language scripts from the choices selected in the menus. These scripts are automatically 
saved as text files named script.tmp , script.log , and error.log located in the directory 
specified by the FL_TMP environment variable. Enter the UNIX™ command “env” in the 
shell window to display the current directory specifications. The script.tmp file is a record 
of the most recent analysis conducted, while script.log is a running record of all analysis 
conducted during the present session. The error.log file is self explanatory. When Flightlab 
is exited, these files are purged and will be started anew the next time Model Analysis is 
utilized. Either of the files can be saved under a different directory and name to be executed 
at a later time in the Model Analysis command line interface or in Gscope. 

The choices for analyzing a model are numerous and not every scenario can be 
covered here. Two examples are provided to illustrate the use of the Model Analysis 
program. The first is a parameter sweep to get cockpit control positions for airspeeds from 
0 to 150 knots in 30 knot increments. The second is a linearization, quasistatic reduction, 
and eigenanalysis of the model using the states [u; w; q; 0; v; p; <J>; r], 
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All choices are made by entering the option number, then pressing “enter”, or placing 
the cursor next to the option with the up and down arrow keys and then pressing “enter”. In 
the following examples, menu titles are displayed in bold capital letters, menu items are 
italicized for ease of understanding. If menus have more than one page, use the “f ’ key to 
scroll to next page and “b” to scroll to previous page. 

• Beginning at MAIN MENU 

• Item 1 Load a Model, select local or central database 

• Select desired model, when finished loading press “e” to return to MAIN MENU 

• Item 4 Test Operations 

• Item 2 Analysis 

• Item 8 Parameter Sweep 

• Item 1 Select Parameter 

• Item 1 Equivalent Airspeed 

• Starting Value 0: Ending Value 150: Incremental Change 30 

• Enter “e” to return to PARAMETER SWEEP OPTIONS 

• Item 2 Select Analysis 

• Item 1 Trim 

• Press “enter” twice, a “y” will be displayed to the right of Trim 

• Enter “e” to return to PARAMETER SWEEP OPTIONS 

• Item 3 Start Sweep, upon completion of calculations, enter “e” to return to 
ANALYSIS SELECTION 

• Enter “e” to return to TEST OPERATIONS 

• Item 4 Plotting Options 
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• Item 1 Current Session 

• Item 6 Parameter Sweep Plot 

• Item 2 Trim Variables 

• Select trim variable, a plot vs airspeed will be displayed 

To obtain a linearized eight state model from the non-linear model with the states [u; w; q; 
0; v; p; (J); r], 

• Beginning at MAIN MENU 

• Item 1 Load a Model, select local or central database 

• Select desired model, when finished loading press “e” to return to MAIN MENU 

• Item 4 Test Operations 

• Item 1 Initialization 

• Item 1 Configuration 

• Item 2 Airframe 

• Item 1 Equivalent vel in I-frame (kts), enter desired velocity at prompt 

• Enter “e”, then 1 to accept choice(s), to return to SELECT A SUBSYSTEM 

• Enter “e” to return to INITIALIZATION 

• Enter “e” to return to TEST OPERATIONS 

• Item 2 Analysis 

• Item 1 Trim, then select 5 to run (aircraft will at trimmed at initial conditions 
specified above. 

• After calculations are complete. Item 3 Linearization 

• Item 4 Quasistatic Reduction, then select option 2 to edit 
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• Item 2 Airframe 

• Select items 22,23,5,12,23,7,10,9, upon completion enter “e”, then 1 to apply and 
return to SELECT A SUBSYSTEM 

• Enter “e” to return to ANALYSIS SELECTION 

• Item 4 Quasistatic Reduction, then enter 4 to execute 

• Item 5 Eigenanalysis 

• Enter “e” to return to TEST OPERATIONS 

• Item 4 Plotting Options 

• Item 1 Current Session 

• Item 1 Eigenvalue Plot 

• Enter “p” to display eigenvalue plot 
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III. SH-60B MODEL DEVELOPMENT 


Development of the SH-60B model involved using both the Model Editor and 
Gscope. Gscope was used to model the control system while the Model Editor was used to 
model the remaining systems. This approach was used because the standard model structures 
present in Model Editor emulated the SH-60B accurately with the exception of the control 
system. 

Each supercomponent has associated with it four script files that are called by the .def 
file when loading the model. The .def file specifics will be discussed shortly. The four script 
files, by convention, have the suffixes .prolog, .exc, .configure, and . epilog and have as the 
prefix the name closely describing the supercomponent group. If selecting standard model 
structures in Model Editor, the four script files are already in place. If modeling a separate 
supercomponent, as is the case here, the script files must be created in conjunction with the 
model development in Gscope. The script files serve the following purposes. The execute 
file, control.exe, is generated in Gscope from the graphical representation of the control 
system as discussed in Chapter II. It is used to instance and connect components and load 
the data associated with each component. The prolog file, control.prolog, contains the data 
needed by the execute file. The epilog file, control.epilog, is where connections are made 
between the control system model created in Gscope an the other supercomponents in the 
model. Finally, the configure file, control .configure, is used to equivalence variables within 
the control system and also let the pilot controls, e.g. cyclic, collective, and pedals, be 
accessed in the Model Analysis program for aircraft dynamic response from control inputs. 

Comments in any of the script files can be included and are preceded by a comment 
marker, //, or by using the “describe” feature in Flightlab. The describe feature allows the 
comments associated with data fields to be displayed when working in Gscope or Model 
Analysis, with a model loaded, by using the describe command. For example, the 
control.prolog file contains the line, k2f=6076.115/3600 “knots to feet per second 
conversion factor”. After going to the world_model_ control_data subgroup, the command 
describe will display: “k2f—knots to feet per second conversion factor”. 




A. MODEL EDITOR PROCEDURES 


Implementation of the SH-60B in the Model Editor was straight forward and 
followed the procedures outlined in Chapter H Data required for the Rotor2, Airframe, 
FlightControl, and Aerolnterference supercomponents were obtained from Sikorsky Aircraft 
in Reference 7. Aerodynamic data required for the Rotorl, main rotor, supercomponent was 
not supplied as Flightlab comes equipped with the necessary tables for the UH-60, which is 
identical to the SH-60B. All fuselage aerodynamic and control system data tables used in 
model construction are listed in Appendix B. 

The flight control supercomponent was selected and saved with the model regardless 
of the fact that a separate flight control system was modeled. This was required to allow the 
.def file to execute the four script files necessary for the flight controls supercomponent. To 
allow the new flight control system to be called and implemented and not the one originally 
selected in Model Editor, the applpath.scr script file was modified. The applpath.scr file is 
called during execution of the .def file and is used to set a search path. In this instance the 
path was the location of the four script files associated with the new flight control system, 
control.prolog, control.exe, control.configure, and control.epilog. 

Upon completion of entering the required data in Model Editor, it was saved and the 
script file sh60b_rigid.def was generated. A file with the .def extension is by convention a 
file that defines a sequence of instruction and other script files to execute [Ref. 2: p. 33], The 
script file sh60b_rigid.def is included in Appendix C. 

Although by examining the sh60b_rigid.def file it can be determined which super and 
sub components were selected and also the data and data file names associated with them, 
it is not an easy task. Immediately following the sh60b_rigid.def in Appendix C is an outline 
which contains all data entered into the Model Editor. It follows a one to one 
correspondence with the Model Editor tree. 

The completed model can now be loaded by executing the sh60b_rigid.def in the 
Gscope window. The execute option is beneath the “File” pull down menu. It is 
recommended that the model first be loaded in Gscope vice going directly to Model Analysis 
and loading it. If errors are encountered during execution, they will be listed as well as the 
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location in the .def file where execution halted. If loading a model in Model Analysis and 
errors are encountered, execution will stop, but no indications are given as to why. 

B. CONTROL SYSTEM MODELING WITH GSCOPE 

1. Model Hierarchy 

A Flightlab model is built using a hierarchial technique. Simple sub-systems are 
built, tested on their own and then combined to form the complete model [Ref. 3: p. 1-1], 
This method allows for a modular design where sub-systems can be easily replaced with 
variations of that sub-system, e.g., a main rotor system with four blades vice three. 

At the lowest level in the hierarchy is a group. A group is made up of a series of 
components and given a name which describes the sub-system of the aircraft. Groups can 
exist within groups. Obviously, a group which contains another group or groups is one level 
higher in the hierarchy. Many sub-system groups are created and linked to form the 
complete aircraft and collectively called the “model” group. The model group is, in turn, 
contained within the “world” group, which is the top level of the hierarchy. 

The remainder of this chapter describes the construction of the SH-60B control 
system and how it is integrated into the part of the model created under Model Editor. The 
presentation will follow a top down order beginning at the top level, which is the world 
group. This is the way models are presented when first opened in Gscope. To aide in ease 
of reading and understanding, all component names will be italicized, group names will be 
bold, names assigned to components will be CAPITALIZED and script file names will be 
underlined . Flightlab does not have an option to name the components as shown in the 
following figures. The names were inserted using a separate graphics utility. 

2. World Group 

The top level in any model is the world group and for the SH-60B model contains 
only the model group. If modeling an entire aircraft, the world group would also contain 
the atmosystem and aerosystem components [Ref. 2: p. 11]. When selected, the model group 
has one group beneath it, the control group. Beneath it are the Sensors, Stabilator, and 
FlightControls groups. Figure 7 shows the progression of the SH-60B model hierarchy. 
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3. Sensors Group 

The sensors group is used to provide state variable input to the various control 
groups during execution of the model. The Sensors group is shown in Figure 8. At the top 
of Figure 8 are two source-transfer function-sink combinations and a single source 
component. The sink components labeled ALPHA and BETA output the models angle of 
attack and sideslip, respectively. The transfer function components are first order low pass 
filters which introduce a very short time delay into the calculations to prevent “cyclic 



data dependency” errors. The source component UNITY provides a constant source of one 
for input into various parts of the control system. 

The aircraft euler angles and euler angle rates are input into the Sensors group via 
the DOFSENSORA source component, while the aircraft body velocities and accelerations 
are input via the DOFSENSORL source component. Since construction of the control 
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system is in Gscope while the remainder of the model was completed under Model Editor, 
the connections of DOFSENSORA and DOFSENSORL to the airframe group must be 
accomplished in the control.epilog file. 

The source components each are connected to bngain components that act as 
demultiplexers to select one of the degrees of freedom from the source component. Each 
blimit component is connected to a bgain component to convert the quantity to the desired 
units, degrees, degrees per second, or knots. Next, the blimit components PLIMIT, 
QLIMrr, RLIMIT, and VDOTLIMIT set upper and lower limits of angular rates and lateral 
acceleration for input into their respective transfer functions. The transfer function 
components labeled PSENSOR, QSENSOR, RSENSOR, and LATACCEL are all second 
order filters which represent the sensor dynamics plus filters in the automatic flight control 
system, (AFCS), and SAS computers. 

4. Stabilator Group 

The incidence angle of the stabilator, horizontal tail, of the SH-60B is programmable, 
the control system for which is presented in Figure 9. The incidence angle is set by the 
control system so as to maintain a nearly level pitch attitude as a function of forward airspeed 
pitch rate, lateral acceleration and collective position, with forward airspeed having the 



Figure 9. Stabilator Group 
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greatest influence on setting the incidence angle. At low airspeeds, below approximately 30 
knots, the stabilator incidence angle is set at 40 degrees, leading edge up, and decreases with 
increasing airspeed. 

Forward airspeed in knots is input into the btransf component VXLAG in Figure 9. 
Collective position in inches is input into the bgain component XCPC. The XCPC gain 
multiplies collective position in inches by ten to convert the quantity to percent collective 
position. These two quantities are then combined to form a column vector by the bsumj 
component STABSUM. The column vector is in turn input into the bsink component 
STABINPUTSINK. 

The bvgain component STABPOSITMAP is a two dimensional lookup table with 
forward airspeed and percent collective position as the two parameters, which are obtained 
from STABINPUTSINK. The output of STABINPUTSINK is not directly connected to 
STABPOSITMAP, but is made in the control.configure file by equivalencing the input and 
output variables of each component. The command: world_ model_control_stabilator_ 
stabpositmap_var=&world_model_control_ stabi!ator_stabinput_y, sets the output of 
STABINPUT equal to the input of STABPOSITMAP. The bvgain component is generally 
used as a scalar or matrix gain block in the feedback loop of a control system. In this case 
the input would be made by directly connecting the components in the model window. 
When the bvgain component is used as a two dimensional lookup table, its input is obtained 
via its “var” parameter. The VECTOR UNITY block supplies a constant input vector of 
[ 1; 1 ]. The second bvgain component GAINMAP operates in the same fashion except it is 
a one dimensional lookup table with the output of VXLAG as its input parameter. 

Lateral acceleration is input into system via the bgain component VYBDOT. The 
gain has a value of one and was used to help clarify this input. Pitch rate is input via the 
bgain component GAIN 1. The btransf component ACTUATOR is second order and 
represents the dynamics of the actuator which drives the stabilator. The blimit component 
RANGELIMIT sets the upper and lower limits of travel at 40 and -8 degrees respectively. 
The value is then converted to radians and sent to the bsink component STABPOSIT. 
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The output of STABPOSIT is equivalenced to the chinge component EITAIL in 
hstabl and hstab2 groups in the control .epilog file. Two connections must made because 
Flightlab models the stabilator as two independent groups which can operate independent 
of each other. The left and right sides of the stabilator are mechanically connected in the 
SH-60B, and therefore operate as one. 

5. Flight Controls Group 

The FlightControls group consists of the Longitudinal, Lateral, Collective, and 
Directional sub-groups. Figure 10, which model the four control axes. These sub-groups 
all operate similarly so only the Longitudinal sub-group will be explained, which is shown 
in Figure 11. 



a. Longitudinal Group 

Excluding the four sub-groups PitchAnalogSAS, PitchDigitalSAS, 
PitchSASInput, and PBA, this sub-group models the mechanical portion of the longitudinal 
axis. The longitudinal control system starts with two source components XBTRM and XB, 
which represent the longitudinal cyclic stick position in inches. The XBTRM component 
is set to a value determined by the trim control matrix which is calculated during the trim 
routine for the model [Ref. 2: p. 27], The XB component is used when flying or performing 
an analysis on the model, otherwise it is set equal to zero. The two values are then summed 
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and sent through the blimiter component B1SLIMIT which sets the upper and lower limits 
of longitudinal cyclic travel at 10 and 0 inches respectively. Inputs from the stability 
augmentation system, SAS, and pitch bias actuator, PBA, are then added. The bgain 
component KXBB1S converts inches of longitudinal cyclic position to degrees of main rotor 
longitudinal pitch angle, commonly called BIS. 

The collective and directional control axes supply inputs to the longitudinal 
control axis via the bgain components GAIN! and GAIN2 respectively. This is known as 
control mixing. The lateral and directional axes also have control mixing, although only 
from the collective axis. The purpose of control mixing is to decouple the helicopter’s 
response by feed forward of control inputs. For example, tail rotor pedal motion is input into 
the longitudinal axis to counter act the pitching moment from the canted tail rotor. Control 
mixing is implimented directly in the control system via the mechanical linkages. 

The sum of collective and tail rotor control mixing inputs is summed and in 
turn summed at SUMJ4 and again summed at SUMJ5 with the control axis bias. The bias 
for each control axis is the angular position of the control surface with the cockpit cyclic, 
pedals and collective in their reference control positions. The SH-60B reference control 
positions are; full aft longitudinal cyclic, full left lateral cyclic, full down collective, and full 
left pedal. For example, the longitudinal axis has a bias of 9.7754 degrees. When all cockpit 
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controls are placed in their reference position, the longitudinal cyclic pitch is 9.7754 degrees. 
The output is then passed through the btransf component PRISERVO, which represent the 
dynamics of the hydraulically actuated primary servos in the helicopter. 

The next portion of the control system sums the output of the PRISERVO 
with the source component B1SCHK. This component is used during development of the 
control system and for inputs directly into the longitudinal pitch angle of the main rotor. The 
control position is converted to radians by the DEG2RAD component and finally input to 
the sinfcBlSOUT. The output of BISOUT is equivalenced to “bis” in the control.configure 
file. This allows the commanded longitudinal cyclic pitch angle to be input into the rotorl 
group. 

b. Stability Augmentation System 

The lateral, longitudinal and directional control axes each have digital and 
analog SAS channels. The collective SAS channel contains the barometric altitude hold, 
radar altimeter hold, and coupler which are autopilot functions, and are not relevant to flight 
dynamic modeling [Ref. 8: Appendix A p. 5]. The analog and digital SAS implemented 
here are not the complete systems present in the SH-60B, but only approximations. The full 
model is proprietary property of Sikorsky Aircraft and therefore could not be supplied for 
this thesis work. 



The pitch analog SAS group in Figure 12 shows the basic pitch SAS 
architecture. Pitch rate is input to the channel via the btransf component PAWASHOUT. 
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The pitch rate signal is then sent through a proportional gain path, PAGAIN2, and a lagged 
path, PALAG. The proportional path provides damping, while the lagged path adds stability 
[ Ref 8: Appendix A p. 5], The five components PAGROUND, ASAS, DSAS, PASUM2, 
and PAGAIN2 make up a logic circuit that determines the gain applied to the signal 
depending on whether the analog and digital SAS channels are on or off. If the digital SAS 
channel is off, the gain applied to the analog SAS signal is 2, and vice versa. This allows an 
increase in control response to help offset the one system being disengaged. Although the 
gain is doubled, the range limits remain the same. If both systems are on the gain is 1, and 
if both are off the gain is 0. The PALIMIT component limits the maximum authority of the 
pitch analog SAS to +/- 6.0 percent, the output of which is summed with the output of the 
pitch digital SAS, refer to Figure 11. 

The pitch digital SAS channel is very close in construction to the analog, so 
reviewing it here would be of little value. Instead, the yaw digital SAS channel will be 
presented. Figure 13 shows the basic architecture of the yaw SAS channel. Yaw rate is fed 
into the yaw SAS channel via the YDWASHOUT component and then sent through a 
proportional and lagged path, for the same reasons as the pitch analog SAS channel. 
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Figure 13. Yaw Digital SAS Group 
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The major difference between the yaw SAS channels, both the digital and 
analog, and the other channels, is a velocity dependent switch which causes a change in 
control laws at 50 knots. Below 50 knots, the yaw rate is fed to proportional and lagged 
channels to enhance yaw damping and stability. Above 50 knots, the yaw channel is used 
to provide turn coordination. In this case, yaw rate, roll rate, and lateral accleration are used 
as feedback to improve the flying qualities. 

Lateral acceleration is fed in through the bgain component YDGAIN3 then 
to a lag filter, YDLAG2. Roll rate is input via the YDGAIN4 component and then summed 
with the lateral acceleration signal at YDSUM4. The bdisw component VX50KTSWITCH1 
serves as the velocity dependent switch. When body velocity in the x-direcition, i.e., forward 
velocity, is greater than 50 knots the switch is set so as to allow the signal to pass. If less 
than 50 knots the switch is set to pass a value of zero, which is input via the bsource 
component ZERO. The group of components at the bottom of Figure 13 serve as a logic 
circuit to set the velocity switch. A value of one is multiplied by a gain of 49 in the 
YDGAIN component and then subtracted from the x-direction body velocity, input via the 
VXKNOTS component, at YDSUM6. The quantity is then input to the YAWSWITCH 
component. If this quantity is less than 1, i.e., velocity less than 50 knots, the 
VX50KTS WITCH 1 is set to pass a signal of zero, if greater than 1, it passes the summed 
lateral acceleration and roll rate signal. The connection between the YAWSWITCH and 
VX50KTSWITCH1 components is accomplished in the control.configure file. The group 
of components in the upper right of Figure 13 serve as a logic circuit to set the applied gain 
to the yaw SAS signal, just as in the pitch analog SAS discussed above. 

The actual tum coordination logic and control law in the yaw axis is 
considerably more complex than modeled here. The full system is considered proprietary 
by Sikorsky Aircraft, and the system modeled here was provided as a functional equivalent 
to the actual system. 



c. Pitch Bias Actuator Group 

The pitch bias actuator, PBA, shown in Figure 14 is a second pitch stability 
channel in addition to the analog and digital SAS channels. The PBA is installed in the SH- 
60B to provide a positive stick position gradient with airspeed. In the bucket airspeed region, 
where longitudinal stick gradient is shallow or slightly reveresed, the PBA applies aft 
longitudinal cyclic requiring forward stick to compensate and creating a positive stick 
position gradient. In addition, pitch rate and pitch attitude are used as feedback to enhance 
the longitudinal damping and stability respectively. 

Pitch rate and theta are input via the blimiter components PBAQUMIT and 
PBATHETLIMIT respectively. Each is multiplied by a gain and then summed at SUMJ1. 
Forward velocity is input via the bvgain component PBA3MP which is a one dimensional 
lookup table that multiplies the velocity by a gain. This output is then summed with the pitch 
rate and theta quantities and again summed with a bias of -4.9 percent. 

The portion of the diagram between SUMJ4 and PBAGAIN3 represents the 
dynamics of the actuator. The first limiter component limits the rate of the actuator to 
7.30 percent per second. The second limiter limits actuator authority to 7.15 percent. The 
PBASWITCH component allows the PBA to be disconnected from the model in order to 
allow an open loop analysis. The PBA cannot be secured in the SH-60B as shown here. The 
last two components combine to convert percent of longitudinal cyclic stick position to 
inches. Refer to Figure 11 for the location of the pitch bias actuator input into the 
longitudinal control axis. 
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IV. SH-60B MODEL ANALYSIS 


A. COMPARISON WITH GENHEL RESULTS 

All of the analyses conducted on the SH-60B for comparison with GenHel utilized 
two script files, trimsweep.scr and powerreg.scr and were done at a gross weight of 19462 
pounds. The trimsweep.scr file outputs cockpit control position, pitch and roll attitude, and 
stabilator incidence angle versus airspeed. The powerreg.scr file gives main rotor and total 
power required for a user specified weight and airspeed range. The outputs of trimsweep.scr 
can all be obtained via the parameter sweep menu in Model Analysis with the exception of 
stabilator incidence angle. Main rotor and total power required versus airspeed cannot be 
obtained in Model Analysis. Both script files are included in Appendix D. 

The first comparison between Flightlab and GenHel results was a level flight trim 
sweep of cockpit control positions versus airspeed. Figure 15. The y-axis limits on each plot 
are full travel of each control axis. 



Longitudinal, lateral and pedal control positions all show good comparison, but the 
collective differs by a nearly constant bias of one inch. This is ten percent, relative to full 
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travel, and was initially deemed unacceptable until a comparison of main rotor power versus 
airspeed was conducted. Figure 16. 



The comparison between Flightlab and GenHel main rotor power required was very 
good. From 0 to 80 knots the difference is negligible, above 80 knots the difference 
increases slightly to 6 percent at 150 knots. Referring back to Figure 15, collective position 
versus airspeed, if the Flightlab curve was moved up so as to be coincident with the GenHel 
curve at 60 knots, the trends between collective position and power required would be 
identical. At low airspeeds Flightlab is slightly higher, both converging to the same value 
at 60 knots and then the difference increasing linearly up to 150 knots. 

Although the collective position differed, the more significant parameter of main 
rotor power agreed well. Flightlab tail rotor plus main rotor power, i.e. total power, is also 
included Figure 16. It shows the expected trend of tail rotor power required decreasing with 
increased airspeed as the vertical tail becomes more effective at providing anti-torque to the 
main rotor. 

The third parameter to compare Flightlab and GenHel models was stabilator 
incidence versus airspeed. Figure 17. Up to 30 knots the results are identical due to the 





Figure 17. Stabilator Incidence Angle vs Airspeed 


that the maximum travel limit of the stabilator is set to 40 degrees. The difference between 
the two increases and then converges to the same value at 150 knots. Part of the difference 
is attributed to the fact that one of the inputs to the stabilator control system is collective 
position. As discussed above, the collective position versus airspeed differed by a nearly 
constant bias of one inch. A second run of stabilator position versus airspeed was conducted 
with a bias of one inch added to the collective input into the stabilator control system. The 
results are shown in Table V, along with original values obtained and the GenHel Results. 

The results with the bias are very close to the GenHel results. The difference between 
collective position for the two programs cannot be explained as GenHel was not available 
for use. It is clear though that the difference between the two are due to computational 
differences within each program as the output variables, power required and stabilator 
incidence angle differed by only a small amount when taking into account the collective bias. 
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EAS (kts) 

GenHel 

Flightlab w/o bias 

Flightlab w/ bias 

0 

40.0 

40.0 

40.0 

30 

40.0 

40.0 

40.0 

60 

20.3 

18.58 

20.13 

90 

4.7 

2.76 

4.33 

120 

3.7 

0.98 

2.58 

150 

0.5 

0.64 

0.64 


Table V. Stabilator Incidence Angle In Degrees vs Airspeed 


The fourth comparison between Flightlab and GenHel was pitch and roll attitude 
versus airspeed. Figure 18. 



Flightlab results showed a slightly higher pitch attitude, theta, over the airspeed range 
than GenHel. Again, part of the difference can be attributed to the collective bias. The 
stabilator incidence angle, which is a factor in pitch attitude, was improved with the 
collective bias accounted for, and so was pitch attitude, but to a smaller extent. Table VI 
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summarizes the results. Roll attitude, phi, compared well up to 30 knots, but diverged from 
the GenHel results. From 60 to 150 knots GenHel showed a level attitude, where Flightlab 
had a left roll angle. 


EAS (kts) 

GenHel 

Flightlab w/o bias 

Flightlab w/ bias 

0 

3.10 

3.48 

3.48 

30 

1.42 

2.28 

2.28 

60 

-0.30 

1.07 

0.79 

90 

-0.71 

0.30 

-.11 

120 

-4.98 

-3.27 

-3.84 

150 

-7.98 

-7.45 

-7.46 


Table VI. Pitch Attitude In Degrees vs Airspeed 


B. CLOSED LOOP ANALYSIS 

The control systems of the model are comprised of the analog and digital SAS 
channels and the pitch bias actuator. State variable feedback is used in these systems to 
provide a stable model, which otherwise is unstable. Both characteristics will be shown in 
this section and the next. All of the analyses therefore conducted in this section are with 
the three systems on, making the model closed loop. All analyses were conducted at a gross 
weight of 19462 pounds, 90 knots equivalent airspeed and at standard sea level conditions. 
The coordinate system for the model is defined with the positive x-axis out the nose of the 
aircraft, positive y-axis out the right side and positive z-axis down. Angular displacements 
are defined in accordance with the right hand rule. 

In Flightlab the nonlinear model can be linearized and expressed in the standard form 
x'=Ax + Bu (1) 

y =Cx + Du (2) 

The A notation is used to denote the closed loop plant. Since state variable feedback is 
implemented in the model; 

A=A-BK (3) 
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where A and B are the open loop plant matrices and K is the gain Matrix. The stability of 
the model can then be determined by finding the eigenvalues of the A matrix. 

The Model Analysis program within Flightlab was used to linearize the model. Prior 
to executing the linearization though, the desired conditions must be set and the model 
trimmed. The result yielded a A matrix that was 77x77, i.e., 77 states, and a B matrix that 
was 77x4. The columns of the B matrix correspond to the four cockpit inputs, longitudinal 
and lateral cyclic, collective and tail rotor pedals. The eigenvalues of the model were all in 
the left hand plane indicating that it was stable. To confirm this, standard control input tests 
were conducted for each of the control axes and the time history response of the aircraft 
states recorded. The nonlinear and linear response are presented for each test. The linear 
responses shown for each of the tests are that of the full 77 state model. 

The first test was a longitudinal cyclic impulse of one inch lasting for three seconds 
with a one second rise and fall time. The time response is given in Figure 19. The linear 



Figure 19. Longitudinal Cyclic Impulse Response 


and nonlinear response correspond very well and dampen out to a steady state value in both 
pitch and forward velocity by approximately 90 seconds. 
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The second test of aircraft stability was a lateral impulse of the same description as 
for the longitudinal test. Refer to Figure 20 for the time history response. The linear 
response has a smaller amplitude in roll and roll rate and dampens out quicker than the 
nonlinear response. The response of the aircraft pitch angle, theta, is also included in the 
figure to demonstrate the roll to pitch coupling of the SH-60B. As a right roll is introduced 
the aircraft nose pitches up and when the lateral cyclic is returned to the trimmed position 
the nose pitches down. The coupling is so strong that the initial response in pitch is nearly 
as high as the roll angle. 



Figure 20. Lateral Cyclic Impulse Response 


To test the directional stability of the model, a tail rotor pedal doublet of one inch 
lasting for three seconds with rise and fall time of one second and a zero delay between each 
impulse was used. Figure 21 shows the time response of the aircraft to such an input. Again 
the linear and nonlinear responses are very close with the linear response having a smaller 
amplitude and damping out quicker than the nonlinear response. 
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Figure 21. Pedal Doublet Response 

The final test of the closed loop response was a collective doublet with the same 
dimensions as the pedal doublet. The response of yaw angle, yaw rate and vertical body 
velocity are shown in Figure 22. In addition to demonstrating the stability in this axis it 
shows the reaction of the aircraft to an increase and then decrease in collective setting. The 
coordinate system for the model has the negative z direction defined as up, therefore the plot 
of Vz body shows the aircraft initially climbing. The corresponding increase in torque 
produced by the main rotor with an increased collective setting causes a yaw to the right, 
defined as positive. Figure 22 actually shows the yaw and yaw rate initially going in the 
opposite direction as would be expected. This is not a discrepancy in the results, but a trait 
of a non-minimum phase system. If at least one pole or zero of the transfer function is in the 
right half plan the system is called non-minimum phase. As none of the poles were in the 
right half plane the transfer function of collective to yaw and yaw rate must have at least one 
zero in the right hand plane. Importing the matrices into MATLAB™ and converting them 
to transfer functions of collective to yaw and yaw rate, it was discovered that each transfer 
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function had 


zero each in the right half plane, thus explaining the yaw and 


response. 



The transfer function of collective to Vz body also had a zero in the right half plane, 
but its effects on the response of Vz body is negligible. This is because the zero was at 118.3 
where for yaw rate it was at 4.65x10 6 . The same can also be said for psi, yaw angle, which 
had a zero at 1.60xl0 4 . 

A linear frequency response analysis was also conducted on the closed loop system 
in the longitudinal and lateral axes. Flightlab can also do a nonlinear frequency response 
analysis, but the differences between the linear and nonlinear time history plots were so 
small, it was assessed that the linear frequency response would also give accurate results. 

The frequency response analysis can be conducted from the Model Analysis program, 
but it does not give the user the option to select the frequency range of interest. The 
frequency range is hard coded in with a range from 1.6 to 100 rad/sec. A second drawback 
is that the preset frequency range is not logarithmically spaced. In order to remedy these two 

43 






drawbacks, the file udf.fun was copied into the users home directory and modified. The file 
is located in the directory $FL_DIR/flightscope/analysis. The command: 
ww=logspace(w 1 ,w2,num), was entered after the original assignment command. The syntax 
for the logspace command: wl is the lowest and w2 is the highest frequency in the range and 
num is the number of elements. The function “Linear FrequencyRespons”' is included in 
Appendix D. Before executing the frequency response analysis in Model Analysis, but after 
loading the model, the command: exec("udf.fun",l), must be executed in the command line 
interface. This allows the udf.fun file in the users home directory to be used in the analysis, 
vice the original one. 

Figure 23 is a frequency response plot of Vx body and pitch angle, theta, for a 
longitudinal cyclic input. At the low frequency spectrum the magnitude of theta is below 



Figure 23. Frequency Response From Longitudinal Cyclic 
Input 


0 db, while that of Vx body starts at 100 db and shows a large amount of damping as 
evidenced by the 80 db/decade decrease in gain. The magnitude crosses the 0 db line at 
approximately 1.5 rad/sec. The interpretation of this plot, as it applies to the physical 




description of the helicopter motion, is that large changes in forward velocity with little 
change in pitch angle are experienced below inputs of 1.5 rad/sec. Above 1.5 rad/sec both 
forward velocity and pitch angle response are highly damped. This type of response is not 
observed in fixed wing aircraft, i.e. the magnitude of the pitch response is higher at lower 
frequencies. Part of the reason is due to the fact that longitudinal, and lateral, inputs are 
made into the main rotor, vice directly into the aircraft body as in the case for a fixed wing 
airplane. It is obvious that an input into a rigid body will have a greater response magnitude 
vice an input into a dynamic system, i.e., main rotor, coupled to the rigid body. The 
magnitude of the airspeed response remains high at low frequencies because the main rotor 
forward and aft pitch, which does have a high gain at low frequencies, acts as a “brake” to 
affect airspeed with little damping. 

Figure 24 is a frequency response plot for Vy body and roll angle, phi, versus lateral 
cyclic. The exact same trend is noted as in the longitudinal case which should be expected. 
Figure 25 is a frequency response plot of pitch angle versus a lateral cyclic input. Notice that 
the magnitude of theta is very close to the magnitude of roll angle for a lateral cyclic input 
due to the cross coupling of the two. This aspect was also discussed in the time history 
response to a lateral input. 






















C. OPEN LOOP ANALYSIS 

The model is open loop with the analog and digital SAS and pitch bias actuator 
disengaged. All three systems can be turned off by setting their flags from one to zero in the 
“initialization/configure/flight controls” menu within Model Analysis. After loading it and 
setting the initial conditions to 19462 pounds gross weight and an equivalent airspeed of 90 
knots the model was linearized. The result was a 77 state plant, the same as with the closed 
loop, but it had three unstable eigenvalues. This 77 state linear model was then reduced to 
a model with six degrees of freedom and eight state variables, [u; w; q; 0; v; p; 4>; r]. The 
reduced order plant also had three unstable eigenvalues, almost exactly in the same location 
as the full 77 state linear model. Figure 26 is an eigenvalue plot of the open loop reduced 
order model. The A and B matrices and the eigenvector and eigenvalues for the reduced 
order model are given in Appendix E for airspeeds of 0, 30, 60, 90, and 120 knots. 
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Figure 26 shows a pair of complex unstable roots and one real unstable root. It will 
be shown that the real root is associated with the longitudinal axis, and that the complex 
roots with the lateral axis. 

In the comparison between the nonlinear and linear response in the previous section, 
the full 77 state plant was used. This section will use the reduced order plant for the control 
input tests for two reasons. It has already been shown that the full linear model gives good 
results when compared to the nonlinear response, so it will be now assessed how the reduced 
order compares with the nonlinear. Also, to determine which modes are associated with the 
unstable eigenvalues in the eight state model. 

The first test was a longitudinal impulse of one inch for a duration of .025 seconds. 
A short duration input was used to simply perturb the model about the trimmed condition and 
observe the response. The closed loop analysis required inputs of a longer duration due to 
the robustness of the stability system. Short impulses, as used in the open loop analysis, had 
no appreciable effect on the aircraft states with the stability system engaged. Figure 27 is the 
time response of pitch angle, theta, and pitch rate to an impulse input. 



Figure 27. Open Loop Response From Longitudinal Impulse 





The linear and nonlinear response are very close up to 15 seconds, after which they 
radically depart. The linear response diverges without oscillation, so it must be associated 
with the positive real eigenvalue. 

The second test conducted was a lateral cyclic impulse also of one inch for a duration 
of .025 seconds. Figure 28 is the time history plot of roll angle, phi, and roll rate. The linear 
and nonlinear responses are very close up to about 20 seconds. After that the two depart, 
but show the same general trend up to 30 seconds. The linear response diverges and 
oscillates, so the lateral mode must be associated with the complex real eigenvalues. In both 
instances, the rapid departure of the linear and nonlinear responses was due to the aircraft 
states departing from the small linear range in which the linear model is valid. 



Figure 28. Open Loop Response From Lateral Impulse 


Control input tests were also conducted for the collective and tail rotor impulse 
inputs, but the responses showed the same instability. This is expected as the strong cross 
coupling in the model, and in helicopters in general, causes a change in the lateral axis when 
an input is made into the collective, for instance. As already demonstrated, even the smallest 
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perturbation caused the model to rapidly depart from the trimmed steady state flight 
condition. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

This thesis presents an outline on the use of Flightlab for rotorcraft modeling and 
analysis. A tutorial was presented to highlight the three programs within Flightlab, Model 
Editor, Gscope, and Model Analysis. The steps to implement a supercomponent, i.e. control 
system, into a model created under Model Editor were also covered. The same generic 
process can be followed for the insertion of any supercomponent that differs from the 
generalized model structures available in Model Editor. 

The results of the SH-60B analysis compared favorably with the results of GenHel, 
Sikorsky Aircrafts modeling and analysis software. All of the comparisons were trimmed 
steady state flight, so the fact that the SAS channels were approximations had little or no 
effect on the results. The largest effect was in the closed loop analysis, but to what extent 
could not be assessed as no similar data was available for the full SAS system responses. 
The closed loop analysis showed that the model was stable as all of the eigenvalues of the 
A matrix were in the left hand plane. It also demonstrated that the linear and nonlinear 
responses of the model to control inputs were very close. The open loop analysis showed 
that the helicopter was unstable in pitch, purely divergent, and in roll, oscillatory and 
divergent. The strong cross coupling between the six degrees of freedom of the linearized 
model was also shown. 

B. RECOMMENDATIONS FOR FURTHER RESEARCH 

The modeling and analysis conducted on the SH-60B yielded good results, but there 
are several improvements that could be made to enhance it’s fidelity. First are the analog and 
digital SAS channels. Both systems implemented for this thesis are not the full scale 
systems, but approximations. Discussed in Chapter ID, these systems are proprietary and 
could not be supplied by Sikorsky Aircraft. If though, through an agreement, the full 
systems could be supplied, they should be implemented. Providing the SAS channels are not 
too complex, modeling them in Flightlab would not be too much effort, and this thesis would 
be of great assistance in doing so. 
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Second, the main rotor was modeled as a rigid system. Flightlab has the capability 
to model it as elastic with up to ten mode shapes. The last two enhancements deal with the 
main rotor system aerodynamics. The induced velocity on the main rotor was modeled with 
the uniform inflow option, which is based on simple momentum theory. Either the free 
vortex or prescribed vortex selective modeling components can be used provided that the 
required data can be obtained. Finally, the airloads were modeled as quasi-steady, the 
dynamic stall option can be used, which allows for more accurate simulation results at and 
near stall conditions on the main rotor. 

Another area for possible future work is the implementation of the aerodynamic data 
for the SH-60B in different mission configurations. Sikorsky Aircraft also supplied the data 
for lift, drag, pitching moment, and aerodynamic interference for toipedoes and external fuel 
tanks on the helicopter. The Model Editor does not contain options for these ‘extra’ 
aerodynamic panels and thus would need to be constructed in Gscope. The procedure would 
be similar to that outlined in this thesis for the control system. 

Including the enhancements listed above would allow the model to be used for 
realistic engineering analysis of changes or improvements to the SH-60B. Suppression of 
high vibration levels is one possible area of research, which is an undesired trait of all 
helicopters and one that has received much attention. One technique for the suppression of 
rotorcraft vibrations is higher harmonic control (HHC), which is an active vibration system. 
One method of HHC is achieved by superimposing 4/rev swashplate motion upon basic 
collective and cyclic control inputs [Ref. 9: p. 3]. Reference 9 is a paper which presents the 
basic theory and results of a HHC system installed on a Army OH-6A helicopter. The paper 
also refers the reader to a plethora of other references for HHC theory. The main rotor is 
already modeled to accept swashplate inputs, as demonstrated, and the control system 
components are flexible enough to model the HHC system. Since HHC is designed to 
suppress fuselage vibration levels, the present modeling of the fuselage would need to be 
looked at in regard to accurately predicting vibration levels since it is a rigid model, and not 
elastic, which would be ideal. 
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Finally, the SH-60B model could be modified to enable real time engineering flight 
simulation for use with the pilots workstation. Flightlab requires the use of a multi-processor 
computer because three programs must be run simultaneously, the visual scene program, the 
model program, and the program to operate the pilots workstation [Ref. 2:p. 52]. Currently, 
the necessary hardware for piloted simulations is present in the computer center visualization 
laboratory at the Naval Postgraduate School. 
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APPENDIX A. FLIGHTLAB COMPONENT CLASS ICONS 



Figure 29. Aerodynamic 
Components 



Figure 30. Analytic 
Components 
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Figure 31. Control 
Components 



Figure 32. Kinematic 
Comonents 
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APPENDIX B. SH-60B DATA FILES 


•Filename sh60fuseclal.sav 
•user Wagner 

♦ SEAHAWK FUSELAGE LIFT VS ALFWF LOW ANGLE MAP 

•MAPNAME =LQFMP 

•MAPTYPE =UVRUVR 

•INPUT =ALFWF 

•OUTPUT =LQF 

CLFAL 13 1 0 

-4.73700e+l -3.91100e+l -2.91000e+l -2.16000e+l 
-9.60000e00 0.95000e00 6.10000e00 9.95000e00 

1.42500e+l 1.90000e+l 2.57300e+l 3.17500e+l 

3.78700e+l 


•Filename sh60fuseclah.sav 
•user Wagner 

• SEAHAWK FUSELAGE LIFT VS ALFWF HIGH ANGLE MAP 

•MAPNAME =LQFMP 

•MAPTYPE =UVRUVR 

•INPUT =ALFWF 

•OUTPUT =LQF 

CLFAH 19 1 0 

-2.10000e+l -2.31000e+l -3.21400e+l -4.28600e+l 
-5.04600e+l -5.21600e+l -4.73700e+l -2.91000e+l 
-9.60000e00 6.10000e00 1.42500e+l 2.57300e+l 

3.78700e+l 4.83900e+l 4.97100e+l 4.62100e+l 

4.06600e+l 3.26300e+l 2.20000e+l 


• Filename sh60fusecmal.saV 
•user Wagner 

• SEAHAWK FUSELAGE PITCH MOMENT VS ALFWF LOW ANGLE MAP 

•MAPNAME =MQFMP 

•MAPTYPE =UVRUVR 
•INPUT =ALFWF 

•OUTPUT =MQF 

CMFAL 13 1 0 

-8.17360e+2 -7.20650e+2 -6.23500e+2 -6.09450e+2 
-4.24390e+2 -2.97340e+2 -1.50760e+2 -3.73800e+l 
1.07500e+l 7.52000e+l 1.14110e+2 1.77610e+2 

2.40430e+2 
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♦Filename sh60fusecmah.sav 
♦user Wagner 

♦ SEAHAWK FUSELAGE PITCH MOMENT VS ALFWF HIGH ANGLE MAP 

♦MAPNAME =MQFMP 

♦MAPTYPE =UVRUVR 

♦INPUT =ALFWF 

♦OUTPUT =MQF 

CMFAH 19 1 0 

-1.32620e+2 -3.81080e+2 -5.97630e+2 
-7.45910e+2 -8.31090e+2 -8.80820e+2 
-8.17360e+2 -6.23500e+2 -4.24390e+2 
-1.50760e+2 1.07500e+l 1.14110e+2 

2.40430e+2 2.95910e+2 2.71910e+2 

1.95290e+2 7.24100e+l -9.92100e+l 

-3.07380e+2 


♦Filename sh60fuseclbal.sav 


♦user Wagner 

♦ SEAHAWK FUSELAGE DELTA LIFT VS PSIWF 

♦MAPNAME =DLQFMP 

♦MAPTYPE =BIV 

♦INPUT1 =PSIWF 

♦INPUT2 =ALFWF 

♦OUTPUT =DLQF 

CLFBAL 13 3 0 


-1.75000e+l 
-1.70000e00 
0.70000e00 
-3.60000e00 
2.20000e00 
-l.OOOOOeOO 
-8.20000e00 
9. SOOOOeOO 
0.OOOOOeOO 
6.2000 0e00 


■1.40000e+l 
-2.lOOOOeOO 
0.20000e00 
•3. OOOOOeOO 
2.90000e00 
•1.lOOOOeOO 
-1.19000e+l 
4.60000e00 
1.2 OOOOeOO 
5.90000e00 


-9.40000e00 
0.OOOOOeOO 
-0.20000e00 
-1.30000e00 
-0..80000e00 
-1.lOOOOeOO 
1.25000e+l 
4.80000e00 
3.70000e00 
4.90000e00 


-4.60000e00 
0.20000e00 
1.40000e00 
0. SOOOOeOO 
0.OOOOOeOO 
-4.60000e00 
1.20000e+l 
2.90000e00 
5.20000e00 


62 














# SEAHAWK FUSELAGE DELTA DRAGVS PSIWF HIGH ANGLE MAP 

#MAPNAME =DDQFMP 

ttMAPTYPE =BIVBIV 

4INPUT1 =PSIWF 

#INPUT2 =ALFWF 

^OUTPUT =DDQF 

CDFBH 19 I 0 

1. 71000e+2 1.68000e+2 1.62000e+2 1.4IOOOe+2 

1.12000e+2 7.60000e+l 2.93000e+l 6.40000e00 

4.60000e00 O.OOOOOeOO 3.00000e00 1.53000e+l 

2.97000e+l 7.60000e+l 1.15000e+2 1.40000e+2 

1.58000e+2 1.68000e+2 1.7X000e+2 






# SEAHAWK FUSELAGE DELTA PITCH MOM VS PSIWF 
#MAPNAME =DMQFMP 


#INPUT1 =PSIWF 

#INPUT2 =ALFWF 

#OUTPUT =DMQF 

CMFBAL 13 3 0 

8.20000e+l 5.90000e+l 4.20000e+l 

2.60000e+l 2.00000e+l 1.80000e+l 

0.OOOOOeOO 0.OOOOOeOO -7.00000e00 
-2.00000e+l -4. 75000e+l -8.70000e+l 
-1.16000e+2 -8.00000e+l -7.60000e+l 
-7.00000e+l -6.20000e+l -3.90000e+l 
-5.OOOOOeOO 0.OOOOOeOO -1.60000e+l 
-3.50000e+l -8.00000e+l -1.06500e+2 
-1.23000e+2 -1.40000e+2 -8.00000e+l 
-9.50000e+l -1.OOOOOe+2 -6.60000e+l 
-5.60000e+l -1.50000e+l 7.50000e+l 

1.20000e+2 1.71000e+2 1.86000e+2 

1.97000e+2 2.05000e+2 2.08000e+2 


#Filename sh60fusecybal.sav 


tMAPNAME =YQFMP 

#MAPTYPE =UVRUVR 

#INPUT1 =PSIWF 

#INPUT2 =ALFWF 


ANGLE 


CYFBAL 13 3 0 

-7.20000e+l -5.95000e+l -4.70000e+l 
-3.47500e+l -2.20000e+l -9.75000e00 
2.50000e00 1.30000e+l 2.35000e+l 

3.75000e+l 5.15000e+l 6.50000e+l 

7.85000e+l -7.20000e+l -5.95000e+l 
-4.70000e+l -3.47500e+l -2.20000e+l 
-9.7 SOOOeOO 2.50000e00 1.30000e+l 

2.35000e+l 3.75000e+l 5.15000e+l 

6.50000e+l 7.85000e+l -7.20000e+l 

-5.95000e+l -4.70000e+l -3.47500e+l 
-2.20000e+l -9.75000e00 2.50000e00 

1.30000e+l 2.35000e+l 3.75000e+l 

5.15000e+l 6.50000e+l 7.85000e+l 
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♦Filename sh60fusecybah.sav 
#user Wagner 

#SEAHAWK FUSELAGE SIDE FORCE VS PSIWF HIGH ANGLE MAP 

#MAPNAME =YQFMP 

#MAPTYPE =UVRUVR 

♦INPUT1 = PSIWF 

♦INPUT2 =PSIWF 

♦OUTPUT =YQF 

CYFBAH 19 3 0 


-4.0000e+l 
-9.6000e+l 
-7.2000e+l 
2.5000e00 
7.8500e+l 
9.5000e+l 
3.9500e+l 
-8.65000e+l 
-8.70000e+l 
-2.20000e+l 
5.15000e+l 
9.95000e+l 
5.95000e+l 
-6.60000e+l 
-9.55000e+l 
-4.70000e+l 
2.35000e+l 
9.40000e+l 
8.00000e+l 


-6.60000e+l 
-9.55000e+l 
-4.70000e+l 
2.35000e+l 
9.40000e+l 
8.00000e+l 
-4.0000e+l 
- 9.6000e+l 
-7.2000e+l 
2.5000e00 
7.8500e+l 
9.5000e+l 
3.9500e+l 
-8.65000e+l 
-8.70000e+l 
-2.20000e+l 
5.15000e+l 
9. 95000e+l 
5.95000e+l 


-8.65000e+l 
-8.70000e+l 
-2.20000e+l 
5.15000e+l 
9.95000e+l 
5.95000e+l 
-6.60000e+l 
-9.55000e+l 
-4.70000e+l 
2.35000e+l 
9.40000e+l 
8.00000e+l 
-4.0000e+l 
-9.6000e+l 
-7.2000e+l 
2.5000e00 
7.8500e+l 
9.5000e+l 
3.9500e+l 


♦Filename sh60fusecnbal.sav 


♦user Wagner 

♦SEAHAWK FUSELAGE FOLLING MOMENT VS PSIWF 

♦MAPNAME =RQFMP 

♦MAPTYPE =UVRUVR 

♦INPUT =PSIWF 

♦OUTPUT =RQF 

CNFBAL 13 3 0 


1.25000e+l 
1.80000e+l 
1.25000e+l 
9.80000e+l 
1.16000e+2 
2.00000e+l 
1.10000e+l 
7.20000e+l 
1.11000e+2 
1. 80000e+l 
1.20000e+l 
4.00000e+l 
1.08000e+2 


1.80000e+l 
1.20000e+l 
4.00000e+l 
1.08000e+2 
1,25000e+l 
1.80000e+l 
1.25000e+l 
9.80000e+l 
1.16000e+2 
2.00000e+l 
1. lOOOOe+l 
7.20000e+l 
1.11000e+2 


2.00000e+l 
1.10000e+l 
7.20000e+l 
1.11000e+2 
1.80000e+l 
1.20000e+l 
4.00000e+l 
1.08000e+2 
1.25000e+l 
1.80000e+l 
1.25000e+l 
9.80000e+l 
1.16000e+2 
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#SEAHAWK FUSELAGE ROLLING MOMENT VS PSIWF HIGH ANGLE MAP 


#MAPNAME =RQFMP 

#MAPTYPE =UVRUVR 



-5.90000e+l -5.10000e+l -4.20000e+l 
-3.10000e+i -2.00000e+l -7.OOOOOeOO 
1.25000e+l 2.00000e+l 1.20000e+l 

1.25000e+l 7,20000e+l 1.08000e+2 

1.16000e+2 1.18000e+2 1.19500e+2 

1.20000e+2 1.20500e+2 1.20750e+2 

1.21000e+2 -5.90000e+l -5.10000e+l 
-4.20000e+l -3.10000e+l -2.00000e+l 
-7.OOOOOeOO 1.25000e+l 2.00000e+l 

1.20000e+l 1.25000e+l 7.20000e+l 

1.08000e+2 1. 16000e+2 1.18000e+2 

1.19500e+2 1.20000e+2 1.20500e+2 

1.20750e+2 1.21000e+2 -5.90000e+l 

-5.10000e+l -4.20000e+l -3.10000e+l 
-2.00000e+l -7.OOOOOeOO 1.25000e+l 
2.00000e+l 1.20000e+l 1.25000e+l 

7.20000e+l 1.08000e+2 1.16000e+2 

1.18000e+2 1.19500e+2 1.20000e+2 

1.20500e+2 1.20750e+2 1.21000e+2 

#Filename sh60fusecrbal.sav 
#user Wagner 

#SEAHAWK FUSELAGE YAWING MOMENT VS PSIWF LOW ANGLE MAP 
#MAPNAME =BQFMP 

#MAPTYPE =BIVBIV 

#INPUT1 =PSIWF 


#OUTPUT =NQF 

CRFBAL 13 3 0 

-1.30000e+2 -1.67000e+2 -1.96000e+2 

-2.03000e+2 -1.56000e+2 -7.80000e+l 

3.50000e+l 1.25000e+2 2.00000e+2 

2.35000e+2 2.46000e+2 2.39000e+2 

2.29000e+2 -9.25000e+l -1.27000e+2 

-1.56000e+2 -1.75000e+2 -1.34000e+2 

2.00000e+l 7.30000e+l 1.48000e+2 

2.42000e+2 2.76000e + 2 2.88000et-2 

2.80000e+2 2.74000e+2 -8-OOOOOetl 

-9.50000e+l -1.00000e+2 -6.60000e+l 

-5.60000e+l -1.5000e+l 7.50000e+l 

1.20000e+2 1.71000e+2 1.86000e+2 

1.97000e+2 2.05000et2 2.08000e+2 


67 













MOMENT 


PSIWF 


ANGLE 


♦ user Wagner 

ISEAHAWK FUSELAGE YAWING 


fMAPNAME 
♦MAPTYPE 
*INPUT1 
♦INPUT2 
♦OUTPUT 

:rfbah 19 3 
4.40000e+2 
2.59000e+2 
-1.30000e+2 
3.50000e+l 
2.2 90 00e+2 
8.80000e+l 
-4.20000e+2 
3.32000e+2 
4.00000e+l 
-1.34000e+2 
2.88000e+2 
1.80000e+2 
-2.08000e+2 
3.92000e+2 
1.60000e+2 
-1.OOOOOe+2 
1.71000e+2 
2.45000e+2 
-4.00000e+l 


=NQFMP 
=BIVBIV 
=PSIWF 
=ALFWF 
=NQF 
0 

3.92000e+2 
1.60000e+2 
-1.96000e+2 
2. OOOOOe+2 
2.45000e+2 
-4.00000e+l 
4.40000e+2 
2.59000e+2 
-9.25000e+l 
7.30000e+l 
2.74000e+2 
8.80000e+l 
-4.2 OOOOe+2 
3.32000e+2 
4.00000e+l 
- 5.60000e+l 
1.97000e+2 
1.80000e+2 
-2.08000e+2 


3.32000e+2 
4.00000e+l 
-1.56000e+2 
2.46000e+2 
1.80000e+2 
-2.08000e+2 
3.92000e+2 
1.60000e+2 
-1.56000e+2 
2.42000e+2 
2.45000e+2 
-4.00000e+l 
4.40000e+2 
2.59000e+2 
-8.00000e+l 
7.50000e+l 
2.08000e+2 
8.80000e+l 
-4.20000e+2 
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((Filename sh60fusecmbah.sav 
#user Wagner 

#SEAHAWK FUSELAGE DELTA PITCH MOM VS PSIWF 

#HI ANGLE MAP 

#MAPNAME =DMQFMP 

(tMAPTYPE =BIV 

#INPUT1 = PSIWF 

#INPUT2 =ALFWF 

#OUTPUT =DMQF 

CMFBAH 19 3 0 

8.20000e+l 8.20000e+l 8.20000e+l 

8.2 0000e+l 8.20000e+l 8.20000e+l 

8.20000e+l 4.20000e+l 2.00000e+l 

0.OOOOOeOO -7 . OOOOOeOO -4.75000e+l 
-1.16000e+2 -1.16000e+2 -1.16000e+2 
-1.16000e+2 -1.16000e+2 -1.16000e+2 
-1.16000e+2 -8.00000e+l -8.00000e+l 
-8.00000e+l -8.00000e+l -8.00000e+l 
-8.00000e+l -8.00000e+l -7.00000e+l 
-3.90000e+l 0. OOOOOeOO -3.50000e+l 
-1.06500e+2 -1.40000e+2 -1.40000e+2 
-1.40000e+2 -1.40000e+2 -1.40000e+2 
-1.40000e+2 -1.40000e+2 -8.00000e+l 
-8.00000e+l -8.00000e+l -8.00000e+l 
-8.00000e+l -8.00000e+l -8.00000e+l 
-1.00000e+2 -5.60000e+l 7.50000e+l 

1.71000e+2 1.97000e+2 2.08000e+2 

2.08000e+2 2.08000e+2 2.08000e+2 

2.08000e+2 2.08000e+2 2.08000e+2 

#Filename sh60fusecmbal.sav 
#user Wagner 

# SEAHAWK FUSELAGE DELTA PITCH MOM VS PSIWF 


ttMAPNAME =DMQFMP 



((INPUT2 =ALFWF 

40UTPUT =DMQF 

CMFBAL 13 3 0 


8.20000e+l 5.90000e+l 4.20000e+l 

2.60000e+l 2.00000e+l 1.80000e+l 

0.OOOOOeOO 0.OOOOOeOO -7.OOOOOeOO 
-2.00000e+l -4.75000e+l -8.70000e+l 
-1.16000e+2 -8.00000e+l -7.60000e+l 
-7.00000e+l -6.20000e+l -3.90000e+l 
-5.OOOOOeOO 0.OOOOOeOO -1.60000e+l 
-3.50000e+l -8,00000e+l -1.06500e+2 
-1.23000e+2 -1.40000e+2 -8.00000e+l 
-9.50000e+l -1.OOOOOe+2 -6.60000e+l 
-5.60000e+l -1.50000e+l 7.50000e+l 

1.20000e+2 1.71000e+2 1.86000e+2 

1.97000e+2 2.05000e+2 2.08000e+2 
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#Filename sh60fusecybal.sav 


#user Wagner 

#SEAHAWK FUSELAGE SIDE FORSE VS PSIWF LO ANGLE MAP 

#MAPNAME =YQFMP 

#MAPTYPE =UVRUVR 

#INPUT1 =PSIWF 

#INPUT2 =ALFWF 

#OUTPUT =YQF 

CYFBAL 13 3 0 


-7.20000e+l 
-3.47500e+l 
2.50000e00 
3.75000e+l 
7.85000e+l 
-4.70000e+l 
-9.75000e00 
2.35000e+l 
6.50000e+l 
-5.95000e+l 
-2.20000e+l 
1.30000e+l 
5.15000e+l 


-5.95000e+l 
-2.20000e+l 
1.30000e+l 
5.15000e+l 
-7.20000e+l 
-3.47500e+l 
2.50000e00 
3.75000e+l 
7.85000e+l 
-4.70000e+l 
-9.75000e00 
2.35000e+l 
6.50000e+l 


-4.70000e+l 
-9.75000e00 
2.35000e+l 
6.50000e+l 
-5.95000e+l 
-2.20000e+l 
1.30000e+l 
5.15000e+l 
-7.20000e+l 
-3.47500e+l 
2.50000e00 
3.75000e+l 
7.85000e+l 


#Filename sh60fusecybah.sav 


fuser Wagner 

fSEAHAWK FUSELAGE SIDE FORCE 
f MAPNAME =YQFMP 

fMAPTYPE =UVRUVR 

fINPUTI =PSIWF 

#INPUT2 =PSIWF - 

fOUTPUT =YQF 

CYFBAH 19 3 0 

-4.0000e+l -6.60000e+l -8.65000e+l 

-9.55000e+l -8.70000e+l 
-4.70000e+l -2.20000e+l 
2.35000e+l 5.15000e+l 

9.40000e+l 9.95000e+l 

8.00000e+l 5.95000e+l 


PSIWF HIGH ANGLE MAP 


-9.6000e+l 
-7.2000e+l 
2.5000e00 

9.5000e+l 


3.9500e+l 
-8.65000e+l 
-8.70000e+l 
-2.20000e+l 
5.15000e+l 
9.95000e + l 
5.95000e+l 
-6.60000e+l 
-9.55000e+l 
-4.70000e+l 
2.35000e+l 
9.40000e+l 
8.00000e+l 


-4.0000e+l 
-9.6000e+l 
-7,2000e+l 
2.5000e00 
7.8500e+l 
9.5000e+l 
3.9500e+l 
-8.65000e+l 
-8.70000e+l 
-2,20000e+l 
5.15000e+l 
9.95000e+l 
5.95000e+l 


-6.60000e+l 
-9.55000e+l 
-4.70000e+l 
2.35000e+l 
9.40000e+l 
8.00000e+l 
-4.0000e+l 
-9.6000e+l 
-7.2000e+l 
2.SOOOeOO 
7.8500e+l 
9.5000e+l 
3.9500e+l 
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sh60htailcll.: 


#user Wagner 

#SEAHAWK HORIZONTAL STAB LIFT COEFFICIENT VS ALFPP1 
# (LOW ANGLE MAP -30 30 DELTA 5) 

#MAPNAME =CLP1MR 

#MAPTYPE =UVSUVS 

#INPUT =ALFPP1 

#OUPUT =CLP1 

CLHL 13 1 0 

-7.45000e-l -7.95000e-l -9.50000e-l -1.01300e00 
-7.45000e-l -3.80000e-l O.OOOOOeOO 3.80000e-l 
7.45000e-l 1.01300e00 9.50000e-l 7.95000e-l 

7.45000e-l 


♦Filename sh60htailclh.sav 
♦user Wagner 

♦ SEAHAWK HORIZONTAL STAB LIFT COEFFICIENT VS ALFPP1 

♦ (HIGH ANGLE MAP -90 90 DELTA 10) 

♦MAPNAME =CLP1MR 

♦ MAPTYPE =UVSUVS 

♦INPUT =ALFPP1 

♦OUPUT =CLP1 

CLHH 19 1 0 

O.OOOOOeOO -1.50000e-l -3.00000e-l -4.40000e-l 
-5.60000e-l -6.50000e-l -7.45000e-l -9.50000e-l 
-7.45000e-l O.OOOOOeOO 7.45000e-l 9.50000e-l 

7.45000e-l 6.50000e-l 5.60000e-l 4.40000e-l 

3.00000e-l 1.50000e-l O.OOOOOeOO 

♦Filename sh60vtailcdl.sav 
♦user Wagner 

♦ SEAHAWK VERTICAL STABILZER DRAG COEFF VS ALFPP3 

♦ (LOW ANGLE MAP -30 30 DELTA 5) 

♦MAPNAME =CDP3MP 

♦MAPTYPE =UVRUVR 
♦INPUT =ALFPP3 

♦OUTPUT =CDP3 

CDVL 13 1 0 

3.60000e-l 2.65000e-l 1.74000e-l 1.18000e-l 

6.60000e-2 3.30000e-2 1.80000e-2 2.10000e-2 

4.40000e-2 9.20000e-2 1.62000e-l 2.48000e-l 

3.55000e-l 
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0.OOOOOeOO -1.50000e-l -3.00000e-l -4.40000e-l 
-5.60000e-l -6.50000e-l -7.45000e-l -9.50000e-l 
-7.45000e-l 0.OOOOOeOO 7.45000e-l 9.50000e-l 

7.45000e-l 6.50000e-l 5.60000e-l 4.40000e-l 

3.OOOOOe-l 1.50000e-l 0.OOOOOeOO 


•Filename sh60vtailcdl.sav 
•user Wagner 

• SEAHAWK VERTICAL STABILZER DRAG COEFF VS ALFPP3 

# (LOW ANGLE MAP -30 30 DELTA 5) 

•MAPNAME =CDP3MP 

•MAPTYPE =UVRUVR 

•INPUT =ALFPP3 

•OUTPUT =CDP3 

CDVL 13 1 0 

3.60000e-l 2.65000e-l 1.74000e-l 1.18000e-l 

6.60000e-2 3.30000e-2 1.80000e-2 2.10000e-2 

4.40000e-2 9.20000e-2 1.62000e-l 2.48000e-l 

3.55000e-l 


•Filename sh60vtailcdh.sav 


• SEAHAWK VERTICAL STABILZER DRAG COEFF VS ALFPP3 

# (HIGH ANGLE MAP -90 90 DELTA 10) 

•MAPNAME =CDP3MP 

•MAPTYPE =UVRUVR 
•INPUT =ALFPP3 

•OUTPUT =CDP3 

CDVH 19 1 0 

1.lOOOOeOO 1.02500e00 9.65000e-l 8.75000e-l 

7.45000e-l 5.75000e-l 3.60000e-l 1.74000e-l 

6.60000e-2 1.80000e-2 4.40000e-2 1.62000e-l 

3.55000e-l 5.80000e-l 7.50000e-l 8.75000e-l 

9.65000e-l 1.02000e00 1.08000e00 
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♦Filename sh60vtailcll.sav 
#user Wagner 

♦ SEAHAWK VERTICAL STABILZER LIFT COEFF VS ALFPP3 

# (LOW ANGLE MAP -30 30 DELTA 5) 

#MAPNAME =CLP3MP 

=UVRUVR 
=ALFPP3 
=CLP3 


-7.30000e-l 
1.60000e-l 
8.90000e-l 


♦MAPTYPE 
#INPUT 
♦OUTPUT 
CLVL 13 1 0 

-l.OOOOOeOO -l.OOOOOeOO -9.30000e-l 
-5.00000e-l -2.80000e-l -6.00000e-2 
3.80000e-l 6.10000e-l 8.20000e-l 

8.90000e-l 


♦Filename sh60vtailclh.sav 
♦user Wagner 

♦SEAHAWK VERTICAL STABILZER LIFT COEFF VS ALFPP3 

# (HIGH ANGLE MAP -90 90 DELTA 10) 

#MAPNAME =CLP3MP 

♦MAPTYPE =UVRUVR 

♦INPUT =ALFPP3 

♦OUTPUT =CLP3 

CLVH 19 1 0 

0.OOOOOeOO -1.20000e-l -2.80000e-l -4.60000e-l 
-6.60000e-l -8.80000e-l -l.OOOOOeOO -9.30000e-l 
-5.00000e-l -6.OOOOOe-2 3.80000e-l 8.20000e-l 

8.900000e-l 8.00000e-l 6.30000e-l 4.80000e-l 

3.20000e-l 1.70000e-l 0.OOOOOeOO 

♦Filename mrint_fusex.sav •• 

♦user Wagner 

♦ SEAHAWK FORE/AFT M.R. DOWNWASH AT FUSELAGE 

♦MAPNAME =EXWFMP 

♦MAPTYPE =BIV 


♦ INPUT2 =AA1FMR 

♦OUTPUT =EKXWF 

RVXFUSEI 11 3 0 

8.OOOOOe-2 1.80000e-l 3.00000e-l 4.30000e-l 

5.50000e-l 6.60000e-l 7.90000e-l 9.00000e-l 

1.03000e00 5.50000e-l 0.OOOOOeOO 0.OOOOOeOO 

1.00000e-l 2.10000e-l 3.20000e-l 4.20000e-l 

5.40000e-l 6.60000e-l 8.00000e-l 9.40000e-l 

5.00000e-l 0.OOOOOeOO -1.20000e-l 2.OOOOOe-2 

8.OOOOOe-2 1.80000e-l 2.80000e-l 4.00000e-l 

' 5.30000e-l 6.70000e-l 8.20000e-l 4.00000e-l 

0.OOOOOeOO 


73 























SAHAWK FORE/AFT M.R. DOWNWASH AT HORIZONTAL TAIL 

I^PNAME = EXPIMP 

^PTYPE =BIV 

NPUT1 =CHIPMR 

STPUT2 =AA1FMR 

LIT PUT =EKXP1 

XHTI 11 3 0 

.OOOOOeOO -2.00000e-l 5.00000e-2 

.00000e-l 5.40000e-l 8.00000e-l 

.04000e00 1.30000e00 1.55000e00 

.00000e-l 0.OOOOOeOO -4.00000e-l 

.00000e-l -2.00000e-l 1.20000e-l 

.60000e-l 6.00000e-l 8.30000e-l 

.06000e00 1.30000e00 6.60000e-l 

.OOOOOeOO -5.60000e-l -8.00000e-l 
.40000e-l -3.20000e-l 4.00000e-2 

.20000e-l 6.00000e-l 8.60000e-l 

.12000e00 5.40000e-l 0.OOOOOeOO 











#Filename mrint_htailz.sav 
#User Wagner 

#SEAHAWK VERTICAL M.R. DOWNWASH AT HORIZONTAL TAIL 

#MAPNAME =EZP1MP 

#MAPTYPE =BIV 

#INPUT1 =CHIPMR 

#INPUT2 =AA1FMR 

#OUTPUT =EKZP1 

# RVZHTI 

RVZHTI 11 3 0 


1.30000e-l 
-1.82000e00 
-1.91000e00 
-1.42000e00 
-9.40000e-l 
-1.98000e00 
-2.14000e00 
-1.35000e00 
-1.91000e00 
-2.14000e00 
-2.16000e00 


-8.00000e-l 
-1.8 6000e00 
-1.94000e00 
-1.14000e00 
-1.84000e00 
-2.04000e00 
-1.89000e00 
-7.80000e-l 
-1.98000e00 
-2.21000e00 
-1.96000e00 


-1.80000e00 
-1.88000e00 
-1.69000e00 
-4.00000e-l 
-1.91000e00 
-2.08000e00 
-1.62000e00 
-1.36000e00 
-2.06000e00 
-2.28000e00 
-1.56000e00 


#Filename fuseint_stabqdyn.sav 
#User Wagner 

#SEAHAWK FDYNAMIC PRESSURE RATIO AT STABILATRO VS ALFWF 
#MAPNAME =QP1MP 

#MAPTYPE =UVR 

#INPUT =ALFWF 

#OUTPUT =QP1QWF 


HTQDYN 1 13 
1.OOOOOeOO 
6.65000e-l 
7.30000e-l 
1.OOOOOeOO 


1.OOOOOeOO 
6.50000e-l 
7.60000e-l 


9.50000e-l 
7.50000e-l 
9.10000e-l 


8.00000e-l 
8.15000e-l 
9.80000e-l 
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#Filename fuseintf_stabvy.sav 
#User Wagner 

#SEAHAWK FUSELAGE LATERAL DOWNWASH ON STABILATOR VS ALFWF 

♦ FVYHT 
FVYHT 110 
0.OOOOOeOO 

♦Filename fuseintf_stabvz.sav 
♦User Wagner 

♦SEAHAWK FUSELAGE DOWNWASH ON STABILATOR VS ALFWF 
♦MAPNAME =EPP1MP 

♦MAPTYPE =UVRUVR 

♦ INPUT = ALFWF 

♦OUTPUT =EPSSP1 

FVZHT 1 13 0 

4.lOOOOeOO 3. SOOOOeOO 2.90000e00 1.90000e00 

-2.50000e-l -1.50000e00 -1.65000e00 -1.40000e00 
-1.95000e00 -4.40000e00 -5.75000e00 -6.65000e00 
-7.50000e00 

♦Filename fuseintf_vertqdyn.sav 
♦User Wagner 

♦ SEAHAWK DYNAMIC PRESSURE RATIO AT VERITCAL TAIL VS PSIWF 

♦MAPNAME =QP3MP 

♦MAPTYPE =BIV 

♦ INPUT1 = PS IWF 


♦OUTPUT =QP3QWF 
VTQDYN 15 3 0 

1.OOOOOeOO 9.75000e-l 9 ; .25000e-l 8.45000e-l 
7.80000e-l 7.20000e-l 6.90000e-l 6.55000e-l 
6.50000e-l 7.60000e-l 7.95000e-l 8.10000e-l 
8.4000 Oe-1 9.20000e-l 1.OOOOOeOO l.OOOOOeOl 
9.55000e-l 8.25000e-l 7.30000e-l 6.60000e-l 
5.90000e-l 5.80000e-l 5.70000e-l 5.80000e-l 
6.10000e-l 6.60000e-l 7.10000e-l 7.80000e-l 
9.10000e-l l.OOOOOeOl 1.OOOOOeOO 9.20000e-l 
7.85000e-l 6.90000e-l 6.15000e-l 5.75000e-l 
5.40000e-l 5.10000e-l 5.20000e-l 5.30000e-l 
5.60000e-l 6.30000e-l 7.30000e-l 8.40000e-l 
1.OOOOOeOO 
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# difference in coordinate system orientation. Flightlab 

# has x forward y out right wing z down where Sikorsky data 

# has x forward y out left wing and z up. 

#MAPNAME =SGP3MP 

#MAPTYPE =BIVBIV 

#INPUT1 =PSIWF 

#INPUT2 =ALFWF 

#OUTPUT =SIGP3 

FVYVT 13 3 0 

-7.80000e00 -7.00000e00 -5.50000e00 -l.OOOOOeOO 
-5.00000e-l 4.00000e-l 1.75000e00 1.90000e00 

2.OOOOOeOO 2.20000e00 4.10000e00 5.30000e00 

5.80000e00 4.35000e00 3.35000e00 2.70000e00 

1.30000e00 1.70000e00 2.80000e00 1.55000e00 

2.00000e-l -1.50000e00 -1.90000e00 -2.65000e00 
-3.55000e00 -4.20000e00 4.40000e00 4.50000e00 

5.OOOOOeOO 5.40000e00 5.20000e00 3.30000e00 

2.20000e00 -9.00000e-l -4.OOOOOeOO -4.10000e00 

-5.20000e00 -6.OOOOOeOO -6.40000e00 
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APPENDIX C. MODEL SCRIPT FILES 


// filename sh60b_rigid.def 

// This file is created by Model Editor 


describe("Model-Id: sh60b_rigid"); 

describe("Model-Description: SH60B Seahawk"); 

exec("$FL_SCRIPTS/startup.scr“,1); 

// First pass instructions: loading data 

goto world; // world 

path("$FL_SCRIPTS/rotorcraft/“); 
pushg (model) 

pushg (data) 
isRotorcraft = 1; 
popg // data 


pushg (data) 

//***•*** begin loading data ******* 

testflag = 0 "flag for test: 0: free fly; 1: wind.tunnel test"; 


















path(”$FL_SCRIPTS/rotorcraft/rotors/blade_element/hub/articulated/"); 
pushg (rotor) 

pushg (data) 
isArticulated = 1; 
popg // data 

pushg (data) 

//***•*** begin loading data ******* 

laghinge = 1 "option for lag hinge: 0: no lag hinge; 1: with lag hinge"; 

flaphof = 1.25 "Rotor flap hinge offset in feet"; 

laghof = 1.25 "Rotor lag hinge offset in feet"; 

feathof = 1.25 "Rotor feathring hinge offset in feet"; 

delta3 = 0 "possitive delta3 for negative pitch-flap coupling"; 

flapkk = 0 "flap hinge rotational spring stiffness"; 

flapcc = 0 "flap hinge damping coefficient"; 

flapangO = 0 "flap hinge spring undeformed angle"; 

//**»**** end loading data ******* 
popg // data 


path("$FL_SCRIPTS/rotorcraft/rotors/blade_element/hub/articulated/lagdam 
per/linear/"); 
pushg (rotor) 

pushg (data) 
isLineardamper = 1; 
popg // data 

pushg (data) 


//**•**** begin loading data ******* 

lagkk = 0 "linear lag damper stiffness"; 

lagcc = 20000 "linear lag damper damping coefficient"; 

lagangO = 0 "linear lag damper spring undeformed angle"; 

//***•*** end loading data ******* 

popg // data 

popg // rotor 


path("$FL_SCRIPTS/rotorcraft/rotors/blade_element/blade/"); 
pushg (rotor) 

pushg (data) 
isBlade = 1; 
popg // data 

pushg (data) 

//******" begin loading data ******* 

firstscpib = 1.25 "Inboard edge location of blade structure, typically 
locates at outmost hinge”; 
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j_SCRIPTS/permpath.s< 


sc (” $FL 




;h("$FL_SCRIPTS/rotorcraft/propu 
;hg (propulsion) 


























































!.Prolog", 1) ; 


ishg (data) 
jpg // data 
jpg // airframe 


»th< n > 

sec(*$FL_SCRIPTS/permpath.scr",1); 

ath("$FL_SCRIPTS/rotorcraft/flight_control/" 
ushg (control) 

ushg (data) 
sFlightControl = 1; 
opg // data 

ushg (data) 

***.*** begin loading data 
******* end loading data ** 
pg // data 


ght_ 
















pushg (fuselage) 


pushg (data) 
isFuseintfVtail = 1; 
popg // data 

pushg (data) 

//******* begin loading data ******* 
vtaostl = [-35 35 5 15]; 

vtaoatl = [-10 10 10 3]; 

load("$FL_LOCAL_DB/tables/rotorcraft//fuseintf_vertqdyn 
vtaost2 = [-90 90 180 -2]; 

vtaoat2 = [-90 90 180 2]; 

load("$FL_LOCAL_DB/tables/rotorcraft//fuseintf_vtailvz. 
vtaost3 = [-30 30 5 13]; 
vtaoat3 = [-10 10 10 3]j 

load("$ FL_LOCAL_DB/tables/rotorcraft//fuseintf_vertvy.s< 
//****•** end loading data ******* 
popg // data 

popg // fuselage 

exec ("applpath.scr", 1) // application path 

pushg (data) 

exec("intffuse.prolog", 1) ; 
popg // data 

popg // fuselage 


exec(”$FL_SCRIPTS/pe 


ith.scr" ,1) ; 
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// file control.exc 

// this file is created in model editor by the model 
// sh60control.mod 

iiintiiiiimniiiiiuiiiiiiiuniimiiinUtiiii't 

GOTO WORLD 
GROUP MODEL 
GROUP CONTROL 
GROUP Sensors 
BLIMITER VDOTLIMIT; 

BLIMITER RLIMIT; 

BLIMITER QLIMIT; 

BLIMITER PLIMIT; 

BSOURCE DOFSENSORA; 

BNGAIN psi; 

BNGAIN theta; 

BNGAIN phi; 

BNGAIN r; 

BNGAIN g; 

BNGAIN p; 

BGAIN psideg; 

BGAIN thetadeg; 

BGAIN phideg; 

BGAIN rdeg; 

BGAIN qdeg; 

BGAIN pdeg; 

BSOURCE ALPHASRC; 

BSOURCE BETASRC; 

BSOURCE UNITY; 

BSOURCE DOFSENSORL; 

BNGAIN U; 

BNGAIN V; 

BNGAIN W; 

BNGAIN VDOT; 

BGAIN UKNOTS; 

BGAIN VKNOTS; 

BGAIN WKNOTS; 

BTRANSF PSENSOR; 

BTRANSF QSENSOR; 

BTRANSF RSENSOR; 

BSINK ALPHA; 

BSINK BETA; 

BTRANSF ALPHADELAY; 

BTRANSF BETADELAY; 

BTRANSF LATACCEL; 

CONNECT RLIMIT(1) TO RSENSOR(1); 

CONNECT VDOTLIMIT(1) TO LATACCEL(1); 

CONNECT VDOT(l) TO VDOTLIMIT(1) ; 

CONNECT DOFSENSORL(1) TO VDOT(1); 

CONNECT QLIMIT(1) TO QSENSOR(1); 

CONNECT qdeg(1) TO QLIMIT(1); 
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CONNECT rdeg(l) TO RLIMIT(l); 

CONNECT PLIMIT(l) TO PSENSOR(l); 
CONNECT pdeg(l) TO PLIMIT(l); 

CONNECT r(1) TO rdeg(l); 

CONNECT q(l) TO qdeg(l); 

CONNECT p(1) TO pdeg(1); 

CONNECT V(l) TO VKNOTS(l); 

CONNECT W(1) TO WKNOTS(l); 

CONNECT DOFSENSORA(1) TO psi(l); 
CONNECT DOFSENSORA(1) TO theta(l); 
CONNECT DOFSENSORA(1) TO phi(l); 
CONNECT DOFSENSORA(1) TO r(1); 
CONNECT DOFSENSORA(1) TO q(1); 

CONNECT DOFSENSORA(1) TO p(l); 
CONNECT psi(l) TO psideg(l); 

CONNECT theta(1) TO thetadeg(l); 
CONNECT phi(1) TO phideg(l); 

CONNECT DOFSENSORL (1) TO U (1 > ; 
CONNECT DOFSENSORL(1) TO V(1) ; 
CONNECT DOFSENSORL (1) TO W(l); 
CONNECT U(1) TO UKNOTS(l); 

CONNECT BETASRC(1) TO BETADELAY(1); 
CONNECT BETADELAY(1) TO BETA(1); 
CONNECT ALPHASRC(1) TO ALPHADELAY(1) ; 
CONNECT ALPHADELAY(1) TO ALPHA(1); 


PARENTG 

GROUP Stabilator 
BGAIN KTS2FPS2; 

BNGAIN DEG2RAD; 

BSINK STABINPUTSINK; 

BNSUMJ STABINPUT; 

BTRANSF VXLAG; 

BVGAIN STABPOSITMAP; 

BSOURCE UNITY; 

BVGAIN GAINMAP; 

BPRODJ PROD1; 

BSUMJ SUMJ1; 

BSINK STABPOSIT; 

BGAIN GAIN1; 

BSUMJ SUMJ2; 

BLIMITER RANGELIMIT; 

BTRANSF ACTUATOR; 

BSOURCE VECTORUNITY; 

BGAIN XCPC; 

CONNECT KTS2FPS2 {1) TO PRODl(l);. 

CONNECT DEG2RAD(1) TO STABPOSIT(1); 
CONNECT RANGELIMIT(1) TO DEG2RAD{1); 
CONNECT STABINPUT(1) TO STABINPUTSINK(1); 
CONNECT VXLAG(1) TO STABINPUT(1); 
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CONNECT STABPOSITMAP(1) TO SUMJ1(2); 
CONNECT GAINMAP(1) TO PRODl(2); 

CONNECT PRODl(l) TO SUMJl(l); 

CONNECT SUMJl(l) TO SUMJ2(1) ; 

CONNECT SUMJ2(1) TO ACTUATOR(1); 

CONNECT GAINl(l) TO SUMJ2(2); 

CONNECT ACTUATOR(1) TO RANGELIMIT(1); 
CONNECT UNITY(1) TO GAINMAP(1); 

CONNECT VECTORUNITY(l) TO STABPOSITMAP (: 
CONNECT XCPC(l) TO STABINPUT(2) ; 


PARENTG 

GROUP FlightControls 
GROUP Collective 
BLIMITER THETLIMIT; 

BSOURCE THOCHK; 

BSUMJ SUMJ4; 

BSINK THETAOOUT; 

BGAIN BIAS; 

BSUMJ SUMJ3 ; 

BGAIN KXCTHO; 

BSOURCE XCILS ; 

BSUMJ SUMJ2; 

BSUMJ SUMJ1; 

BSOURCE XC; 

BSOURCE XCTRM; 

BGAIN DEG2 RAD ; 

BTRANSF PRISERVO; 

CONNECT THOCHK(1) TO SUMJ4(2); 
CONNECT SUMJ4(1) TO DEG2RAD(1); 
CONNECT XC (1) TO SUMJ1 (1) ; 

CONNECT XCTRM(1) TO SUMJ1(2); 
CONNECT SUMJ1 (1) TO THETLIMIT (1) ; 
CONNECT THETLIMIT(1) TO SUMJ2(1); 
CONNECT XCILS (1) TO SUMJ2(2); 
CONNECT SUMJ2(1) TO KXCTHO(1); 
CONNECT KXCTHO(1) TO SUMJ3(1); 
CONNECT BIAS(1) TO SUMJ3(2); 

CONNECT DEG2RAD (1) TO THETAOOUT (1) ; 
CONNECT SUMJ3 (1) TO PRISERVO (1); 
CONNECT PRISERVO(1) TO SUMJ4(1); 


PARENTG 

GROUP Longitudinal 
BSUMJ SUMJ4; 

BSUMJ SUMJ8; 
BSOURCE B1SCHK; 
BSUMJ SUMJ6; 
BSOURCE XB; 


114 



BSUMJ SUMJ1; 

BLIMITER B1SLIMIT; 

BSUMJ SUMJ3; 

BGAXN KXBB1S; 

BSUMJ SUMJ5 ; 

BGAIN BIAS; 

BSINK BlSOUT; 

BGAIN DEG2RAD; 

BSUMJ SUMJ7; 

GROUP PitchAnalogSAS 
BGAIN PAGROUND; 

BDISW DSAS; 

BDISW ASAS; 

BSUMJ PASUM2; 

BGAIN PAGAIN3; 

BLIMITER PALIMIT; 

BTRANSF PALAG; 

BSUMJ PASUMl; 

BGAIN PAGAIN1; 

BGAIN PAGAIN2; 

BTRANSF PAWASHOUT; 

CONNECT PASUM2 (1) TO PALIMIT (1); 
CONNECT ASAS (1) TO PAGAIN3 (1) ; 
CONNECT PAGAIN3 (1} TO PASUM2 (1) ; 
CONNECT PAGROUND (1) TO DSAS (1) ; 
CONNECT PAGROUND(1) TO ASAS(1) ; 
CONNECT ASAS(1) TO DSAS(2); 

CONNECT DSAS(1) TO PASUM2(2>; 
CONNECT PALAG(1) TO PASUMl(2); 
CONNECT PASUMl(1) TO ASAS(2); 
CONNECT PAGAIN1 (1) TO PALAG(1); 
CONNECT PAGAIN2(1) TO PASUMl(1); 
CONNECT PAWASHOUT (1) TO PAGAIN1 (1) ; 
CONNECT PAWASHOUT(1) TO PAGAIN2(1) ; 


PARENTG 

GROUP PitchDigitalSAS 
BGAIN PDGROUND; 

BDISW ASAS; 

BDISW DSAS; 

BSUMJ PDSUM2; 

BGAIN PDGAIN3; 

BLIMITER PDLIMIT; 

BSUMJ PDSUM1; 

BTRANSF PDLAG; 

BGAIN PDGAIN1; 

BGAIN PDGAIN2; 

BTRANSF PDWASHOUT; 

CONNECT PDGAINl(l) TO PDLAG(1); 
CONNECT PDWASHOUT(1) TO PDGAIN1(1); 
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CONNECT PDSUM2(1) TO PDLIMIT(1); 
CONNECT PDGROUND(1) TO ASAS(1); 
CONNECT PDGROUND(1) TO DSAS(1); 
CONNECT DSAS(l) TO ASAS(2); 

CONNECT ASAS(1) TO PDSUM2(2); 
CONNECT DSAS(l) TO PDGAIN3(1); 
CONNECT PDGAIN3(1) TO PDSUM2(1) ; 
CONNECT PDSUMl(l) TO DSAS(2); 
CONNECT PDGAIN2(1) TO PDSUM1(2); 
CONNECT PDLAG(1) TO PDSUM1(1); 
CONNECT PDWASHOUT(1) TO PDGAIN2(1); 


PARENTG 

GROUP PitchSASInput 
BGAIN DEG2IN; 

BGAIN PER2DEG; 

BLIMITER PSASLIMIT; 

BTRANSF PSASACTUATOR; 

CONNECT PSASLIMIT(1) TO PER2DEG(1) ; 
CONNECT PER2DEG(1) TO DEG2IN(1); 

CONNECT PSASACTUATOR(1) TO PSASLIMIT(1); 


PARENTG 
BSUMJ SUMJ2; 

GROUP PBA 
BSOURCE ZERO; 

BDISW PBASWITCH; 

BGAIN PBABIAS; 

BLIMITER PBATHETLIMIT; 

BLIMITER PBAQLIMIT; 

BGAIN PER2DEGS; 

BLIMITER PBALIMIT2; 

BLIMITER PBALIMIT1; 

BINTGR INTEGRATOR; 

BGAIN PBAGAIN3; 

BSUMJ SUMJ4; 

BSUMJ SUMJ3; 

BSUMJ SUMJ1; 

BGAIN PBAGAIN2; 

BGAIN PBAGAIN1; 

BGAIN DEG2IN; 

BSUMJ SUMJ2 ; 

BVGAIN PBA3MP; 

BSOURCE UNITY; 

CONNECT ZERO(1) TO PBASWITCH(1); 
CONNECT PBASWITCH(1) TO PER2DEGS(1); 
CONNECT PBALIMIT2(1) TO PBASWITCH(2); 
CONNECT PER2DEGS (1) TO DEG2INU); 
CONNECT PBAQLIMIT(1) TO PBAGAIN2(1); 
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CONNECT PBATHETLIMIT(1) TO PBAGAINl(l); 
CONNECT SUMJ3(1) TO SUMJ4(1); 

CONNECT PBABIAS(1) TO SUMJ3(2); 

CONNECT SUMJ4(1) TO PBAGAIN3(1); 

CONNECT PBAGAIN3(1) TO PBALIMIT1(1); 
CONNECT PBALIMITl(1) TO INTEGRATOR(1); 
CONNECT INTEGRATOR(1) TO PBALIMIT2(1); 
CONNECT PBALIMIT2(1) TO SUMJ4(2); 
CONNECT PBAGAIN1(1) TO SUMJ1(2); 

CONNECT PBAGAIN2{1) TO SUMJl(l); 

CONNECT SUMJ1(1) TO SUMJ2(1); 

CONNECT SUMJ2(1) TO SUMJ3(1); 

CONNECT PBA3MP(1) TO SUMJ2(2); 

CONNECT UNITY(1) TO PBA3MP(1); 


PARENTG 
BSOURCE XBTRM; 

BGAIN GAIN1; 

BGAIN GAIN2; 

BTRANSF PRISERVO; 

CONNECT SUMJ8(1) TO SUMJ4{2); 

CONNECT SUMJ4(1) TO SUMJ5(1); 

CONNECT KXBB1S(1) TO SUMJ4(1); 

CONNECT SUMJ3(1) TO KXBB1S(1); 

CONNECT GAIN2(1) TO SUMJ8(2); 

CONNECT GAINl(l) TO SUMJ8(1); 

CONNECT B1SCHK(1) TO SUMJ6(2); 

CONNECT SUMJ6(1) TO DEG2RAD(1); 

CONNECT BIAS(1) TO SUMJ5(2); 

CONNECT XBTRM(1) TO SUMJ1(2); 

CONNECT XB(1) TO SUMJ1I1); 

CONNECT DEG2RAD(1) TO BlSOUT(l); 

CONNECT PitchSASInput_DEG2IN(l) TO SUMJ3(2); 
CONNECT PitchAnalogSAS_PALIMIT(1) TO SUMJ7(1); 
CONNECT PitchDigitalSAS_PDLIMIT(l) TO SUMJ7(2); 
CONNECT SUMJl(l) TO B1SLIMIT(1); 

CONNECT B1SLIMIT(1) TO SUMJ2(2); 

CONNECT SUMJ2(1) TO SUMJ3(1); 

CONNECT PBA_DEG2IN(1) TO SUMJ2(1}; 

CONNECT SUMJ5(1) TO PRISERVO(1); 

CONNECT PRISERVO(1) TO SUMJ6(1); 

CONNECT SUMJ7(1) TO PitchSASInput_PSASACTUATOR (1) ; 


PARENTG 
GROUP Lateral 
BSOURCE A1SCHK; 
BSUMJ SUMJ5; 
BSUMJ SUMJ1; 
BSOURCE XATRM; 
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BSOURCE XA; 

BLIMITER A1SLIMIT; 

BSUMJ SUMJ2; 

BGAIN KXAA1S; 

BSUMJ SUMJ4; 

BGAIN BIAS; 

BSINK A1SOUT; 

BGAIN DEG2RAD; 

GROUP ROLLANALOGSAS 
BGAIN RAGROUND; 

BGAIN RAGAIN4; 

BSUMJ RASUM3; 

BDISW ASAS; 

BDISW DSAS; 

BLIMITER RALIMIT3; 

BTRANSF RALAG; 

BSUMJ RASUM1; 

BGAIN RAGAIN3; 

BSUMJ RASUM2; 

BGAIN RAGAIN2 ; 

BLIMITER RALIMIT2; 

BGAIN RAGAIN1; 

BTRANSF RALAG2; 

BLIMITER RALIMIT1; 

CONNECT RASUM3 (1) TO RALIMIT3(1) 
CONNECT DSAS(1) TO RASUM3(2); 
CONNECT ASAS(l) TO DSAS (2); 
CONNECT RAGROUND(1) TO ASAS(l); 
CONNECT RAGROUND! 1) TO DSAS(l); 
CONNECT RAGAIN4(1} TO RASUM3(1); 
CONNECT ASAS(1) TO RAGAIN4(1); 
CONNECT RAGAIN1(1) TO RALAG(1); 
CONNECT RALAG (1) TO RASUM1 (2 ) ; 
CONNECT RASUM1 (1) TO RASUM2 (1) ,- 
CONNECT RASUM2 (1) TO ASAS (2) ; 
CONNECT RALIMIT2(1) TO RAGAIN3(1 
CONNECT RAGAIN3(1) TO RALAG2(1) ; 
CONNECT RALAG2(1) TO RASUM2(2); 
CONNECT RAGAIN2(1) TO RASUMl(l); 
CONNECT RALIMIT1(1) TO RAGAIN2(1 


PARENTG 

GROUP ROLLDIGITALSAS 
BGAIN RDGROUND; 

BGAIN RDGAIN4; 

BSUMJ RDSUM3; 

BDISW DSAS; 

BDISW ASAS; 

BLIMITER RDLIMIT2; 
BSUMJ RDSUM2; 




BGAIN RDGAIN3; 

BSUMJ RDSUM1; 

BTRANSF RALAG; 

BLIMXTER RDLIMITl; 

BGAIN RDGAIN2; 

BGAIN RDGAIN1; 

CONNECT RDSUM3{1) TO RDLIMIT2(1); 
CONNECT RDGAIN4(1) TO RDSUM3(1) ; 
CONNECT DSAS(1) TO RDGAIN4(1); 
CONNECT ASAS(1) TO RDSUM3(2); 
CONNECT DSAS(l) TO ASAS(2); 
CONNECT RDGROUND(1) TO DSAS(1 ); 
CONNECT RDGROUND(1) TO ASAS(1); 
CONNECT RDLIMITl(1) TO RDGAIN3(1); 
CONNECT RDSUM1 (1) TO RDSUM2 (1) ; 
CONNECT RDGAIN3(1) TO RDSUM2(2); 
CONNECT RDGAIN2(1) TO RDSUM1(1); 
CONNECT RALAG(l) TO RDSUMl(2); 
CONNECT RDGAIN1(1) TO RALAG(1); 
CONNECT RDSUM2(1) TO DSAS(2) ; 
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CONNECT SUMJ3(1) TO 
CONNECT GAINl(l) TO 
CONNECT SUMJ 4(1) TO 
CONNECT PRISERVO(1) 
CONNECT SUMJ6(1) TO 


SUMJ4(1); 

SUMJ3(2); 

PRISERVO(1); 

TO SUMJ5(1); 

ROLLSASINPUT_RSASACTUATOR(l); 


PARENTG 

GROUP Directional 
BSUMJ SUMJ6; 

GROUP YAWSASINPUT 
BGAIN DEG2IN; 

BGAIN PER2DEG; 

BLIMITER YAWSASLIMIT; 

BTRANSF YSASACTUATOR; 

CONNECT YAWSASLIMIT{1) TO PER2DEG(1); 
CONNECT PER2DEG(1) TO DEG2IN(1); 

CONNECT YSASACTUATOR (1) TO YAWSASLIMIT (1) ; 


PARENTG 

GROUP YAWDIGITALSAS 
BSUMJ YDSUM4; 

BTRANSF YDWASHOUT; 

BGAIN YDGAIN3; 

BTRANSF YDLAG2; 

BGAIN YDGAIN4; 

BGAIN YDGAIN2; 

BSOURCE ZERO; 

BDISW VX50KTSWITCH1; 

BTRANSF YDLAG1; 

BSUMJ YDSUM1; 

BGAIN YDGAIN1; 

BDISW ASAS; 

BDISW DSAS; 

BSUMJ YDSUM3; 

BGAIN YDGAIN5; 

BGAIN YDGROUND; 

BSUMJ YDSUM2; 

BSOURCE YDUNITY; 

BSINK YAWSWITCH; 

BSUMJ YDSUM6; 

BGAIN YDGAIN; 

BGAIN VXKNOTS; 

BLIMITER YDLIMIT; 

CONNECT VX50KTSWITCH1(1) TO YDSUM2(1); 
CONNECT ZERO(1) TO VX50KTSWITCH1(1) ; 
CONNECT YDSUM4U) TO VX50KTSWITCH1 (2) ; 
CONNECT YDLAG2 < 1) TO YDSUM4 (1) ; 
CONNECT YDGAIN4(1) TO YDSUM4(2); 
CONNECT YDLAG1(1) TO YDSUM1(1); 
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CONNECT YDGAIN2(1) TO YDLAG1(1); 
CONNECT YDWASHOUT(1) TO YDGAIN2(1); 
CONNECT YDWASHOUT(1) TO YDGAIN1(1); 
CONNECT YDGAIN3(1) TO YDLAG2(1) ; 
CONNECT YDSUM1(1) TO YDSUM2(2); 
CONNECT YDGAIN1(1) TO YDSUM1(2) ; 
CONNECT YDGROUND(1) TO ASAS{1); 
CONNECT YDGROUND(1) TO DSAS(1); 
CONNECT DSAS {1) TO ASAS(2); 

CONNECT ASAS(1) TO YDSUM3(2); 
CONNECT DSAS(1) TO YDGAIN5(1); 
CONNECT YDGAIN5 (1) TO YDSUM3 (1) ; 
CONNECT YDSUM2(1) TO DSAS(2); 
CONNECT YDUNITY(1) TO YDGAIN(l); 
CONNECT YDSUM6(1) TO YAWSWITCH(1); 
CONNECT VXKNOTS(l) TO YDSUM6(2); 
CONNECT YDGAIN(l) TO YDSUM6(1); 
CONNECT YDSUM3(1) TO YDLIMIT(1); 


PARENTG 

BSOURCE THRCHK; 

BSUMJ SUMJ5; 

BSUMJ SUMJ1; 

BSOURCE XPTRM; 
BSOURCE XP; 

BLIMITER THETTRLIMIT; 
BSUMJ SUMJ2; 

BGAIN KXPTHR; 

BSUMJ SUMJ4; 

BGAIN BIAS; 

BSINK THETTROUT; 
BGAIN DEG2RAD; 

GROUP YAWANALOGSAS 
BTRANSF YAWASHOUT; 
BSOURCE ZERO; 

BGAIN YAGAIN4; 

BSUMJ YASUM2; 

BDISW VX50KTSWITCH2; 
BGAIN YAGAIN3; 

BDISW VX50KTSWITCH1; 
BLIMITER YALIMIT3; 
BDISW DSAS; 

BDISW ASAS; 

BSUMJ YASUM3; 

BGAIN YAGAIN5; 

BGAIN YAGROUND; 

BSUMJ YASUM1; 

BGAIN YAGAIN2; 

BGAIN YAGAIN1; 

BTRANSF YALAG; 
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BLIMITER YALIMIT1; 

BLIMITER YALIMIT2; 

CONNECT ZERO(1) TO VX50KTSWITCH2(1); 
CONNECT VX50KTSWITCH2(1) TO YASUM2(2); 
CONNECT YASUM2(1) TO ASAS(2); 

CONNECT YASUMl(l) TO YASUM2(1); 

CONNECT YALAG (1) TO YASUM1 (2) ; 

CONNECT YAGAIN1(1) TO YASUMl(l); 

CONNECT YAGAIN4(1) TO VX50KTSWITCH2(2); 
CONNECT VX50KTSWITCH1(1) TO YALAG(1); 
CONNECT YAGAIN2 (1) TO VX50KTSWITCH1 (1) ; 
CONNECT YAGAIN3 (1) TO VX50KTSWITCH1 (2) ; 
CONNECT YALIMIT2(1) TO YAGAIN2{1) ; 
CONNECT YAGROUND(1) TO DSAS(l); 

CONNECT YAGROUND(1) TO ASAS(l); 

CONNECT ASAS(1) TO DSAS(2); 

CONNECT DSAS(1) TO YASUM3(2); 

CONNECT ASAS(1) TO YAGAIN5(1); 

CONNECT YAGAIN5 (1) TO YASUM3(1); 

CONNECT YASUM3(1) TO YALIMIT3(1); 
CONNECT YAWASHOUT(1) TO YALIMIT1(1) ; 
CONNECT YAWASHOUT(1) TO YALIMIT2(1) ; 
CONNECT YALIMIT1(1) TO YAGAIN1(1) ; 


PARENTG 
BSUMJ SUMJ3; 

BGAIN GAIN1; 

BTRANSF PRISERVO; 

CONNECT THRCHK(1) TO SUMJ5(2); 

CONNECT SUMJ5(1) TO DEG2RAD(1); 

CONNECT YAWSASINPUT_DEG2IN(1) TO SUMJ2(2); 
CONNECT YAWANALOGSAS_YALIMIT3 (1) TO SUMJ6(1); 
CONNECT BIAS(1) TO SUMJ4(2) ; 

CONNECT THETTRLIMIT(1) TO SUMJ2(1); 

CONNECT SUMJl(l) TO THETTRLIMIT(1); 

CONNECT XP(1) TO SUMJl(l); 

CONNECT XPTRM(l) TO SUMJ1(2); 

CONNECT DEG2RAD(1) TO THETTROUT(1) ; 

CONNECT SUMJ2(1) TO KXPTHR(1); 

CONNECT KXPTHR(1) TO SUMJ3(1); 

CONNECT SUMJ3(1) TO SUMJ4(1); 

CONNECT GAINl(l) TO SUMJ3(2); 

CONNECT SUMJ4(1) TO PRISERVO(1); 

CONNECT PRISERVO(1) TO SUMJ5(1); 

CONNECT SUMJ6(1) TO YAWSASINPUT_YSASACTUATOR (1) ; 
CONNECT YAWDIGITALSAS_YDLIMIT (1) TO SUMJ6 ( 2 ) ; 


PARENTG 

CONNECT Collective_KXCTHO (1) TO Longitudinal_GAINl ( 
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_RSENSOR (1 


FlightControls_Directional_YAWDIGITALSAS_YDWASHOUT (1) ; 

CONNECT Sensors_RSENSOR (1) TO 

FlightControls_Directional_YAWANALOGSAS_YAWASHOUT(1); 

CONNECT Sensors_LATACCEL(l) TO 

FlightControls_Directional_YAWDIGITALSAS_YDGAIN3(1); 

CONNECT Sensors_PSENSOR(l) TO 

FlightControls_Directional_YAWDIGITALSAS_YDGAIN4(1); 

CONNECT Sensors_PSENSOR(1) TO 

FlightControls_Directional_YAWANALOGSAS_YAGAIN4 (1) ; 

CONNECT Sensor s_LATACCEL (1) TO 

FlightControls_Directional_YAWANALOGSAS_YAGAIN3(1); 

CONNECT Sensors_LATACCEL(1) TO Stabilator_KTS2FPS2(1); 

CONNECT Sensors_UNITY(1) TO 

FlightControls_Lateral_ROLLANALOGSAS_RAGROUND (1) ; 

CONNECT Sensors_UNITY (1) TO 

FlightControls_Lateral_ROLLDIGITALSAS_RDGROUND(1) ; 

CONNECT Sensors_UNITY(1) TO 

FlightControls_Longitudinal_PitchDigitalSAS_PDGROUND(1); 

CONNECT Sensors_UNITY (1) TO 

FlightControls_Longitudinal_PitchAnalogSAS_PAGROUND(l); 

CONNECT Sensors_UKNOTS(1) TO Stabilator_VXLAG(1); 

CONNECT Sensors_QSENSOR (1 }■ TO Stabilator_GAINl (1) 

CONNECT Sensors_UNITY(l) TO FlightControls_Longitudinal_BIAS(1); 
CONNECT Sensors_UNITY(1) TO FlightControls_Lateral_BIAS(1); 
CONNECT Sensors_UNITY(1) TO FlightControls_Collective_BIAS(1); 
CONNECT Sensors_UNITY(1) TO FlightControls_Directional_BIAS(1) ; 
CONNECT Sensors_thetadeg(1) TO 

FlightControls_Longitudinal_PBA_PBATHETLIMIT (1) ; 

CONNECT Sensors_UNITY(l) TO FlightControls_Longitudinal_PBA_PBAB - ' 
CONNECT Sensors_QSENSOR(1) TO 
FlightControls_Longitudinal_PBA_PBAQLIMIT(1); 

CONNECT Sensors_UNITY(l) TO 

FlightControls_Directional_YAWANALOGSAS_YAGROUND(1) ; 

CONNECT Sensors_UNITY(1) TO 

FlightControls_Directional_YAWDIGITALSAS_YDGROUND(l) ; 

CONNECT Sensors_PSENSOR(1) TO 

FlightControls_Lateral_ROLLANALOGSAS_RAGAINl(1); 

CONNECT Sensors_phideg(1) TO 

FlightControls_Lateral_ROLLANALOGSAS_RALIMIT2(1); 

CONNECT Sensors_PSENSOR(1) TO 

FlightControls_Lateral_ROLLDIGITALSAS_RDGAINl(1); 

CONNECT Sensors_PSENSOR(1) TO 

FlightControls_Lateral_ROLLDIGITALSAS_RDGAIN2 (1) ; 





CONNECT 


Sensors_phideg(1) TO 
FlightControls_Lateral_ROLLDIGITALSAS_RDLIMITl(1); 

CONNECT Sensors_QSENSOR(1) TO 

FlightControls_Longitudinal_PitchAnalogSAS_PAWASHOUT (1) ; 
CONNECT Sensor s__QSENSOR(l) TO 

FlightControls_Longitudinal_PitchDigitalSAS_PDWASHOUT(1); 
CONNECT Sensors_PSENSOR(l) TO 

F1ightCont ro1s_Lateral_ROLLANALOGSAS_RALIMITl(1); 

CONNECT Sensors_UKNOTS(1) TO 

FlightControls_Directional_YAWDIGITALSAS_VXKNOTS(1); 


PARENTG 


PARENTG 


GOTO MODEL 
GOTO CONTROL 
GOTO Sensors 

VDOTLIMIT_N =1 
VDOTL IMIT_LL 
VDOTLIMIT_UL 

RLIMIT_N =1 
RLIMIT_LL =-40.0; 

RLIMIT_UL =40.0; 

QLIMIT_N = 1; 

QLIMIT_LL =-40.0; 

QLIMIT_UL =40.0; 

PLIMIT_N =1; 

PLIMIT_LL =-40.0; 

PLIMIT_UL =40.0; 

DOFSENSORA_N =9; 

DOFSENSORA_U =&WORLD_DATA_PFRASENS ; 


psi_NI =9; 



phi_NI =9; 

phi_NO =1; 

phi_K =[100000000]; 


=-16.085; 

=16.085; 
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W_NO Bill 

W_K =[000001000]; 



VDOT_K = [00000001 0]; 

UKNOTS_N = 1; 

UKNOTS_K =&WORLD_MODEL_CONTROL_DATA_F2K; 
VKNOTS_N = 1; 

VKNOTS_K = &WORLD_MODEL_CONTROL_DATA_F 2 K ; 

WKNOTS_N =1; 

WKNOTS_K =6cWORLD_MODEL_CONTROL_DATA_F2K; 

PSENSOR_N =2; 

P S ENSOR__NUM = &WORLD_MODEL_CONTROL_DATA_SNUM ; 
PSENSOR_DEN =&WORLD_MODEL_CONTROL_DATA_SDEN ; 

QSENSOR_N =2 ; 

QSENSOR_NUM =&WORLD_MODEL_CONTROL_DATA_SNUM ; 
Q S ENSOR_DEN = &WORLD_MODEL_CONTROL_DATA_SDEN ; 

RSENSOR_N =2; 

RSENSOR_NUM =&WORLD_MODEL_CONTROL_DATA_SNUM ; 
RSENSOR_DEN =&WORLD_MODEL_CONTROL_DATA_SDEN ; 

ALPHA_N =1; 


ALPHADELAY_N =1; 

ALPHADELAY_NUM =[0 1] ; 

ALPHADELAY_DEN =[.0001 1]; 

BETADELAY_N =1; 

B ETADELAY_NUM = [ 0 1 ] ; 

BETADELAY_DEN =[.0001 1] ; 

LATACCEL_N =2 ; 

LATACCEL_NUM =&WORLD_MODEL_CONTROL_DATA_LATACELNUM ; 

LATACCEL_DEN =&WORLD_MODEL_CONTROL_DATA_LATACELDEN ; 



KTS2FPS2_N =1; 
KTS2FPS2_K =1; 
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RANGEL IMIT_N cl; 

RANGELIMIT_LL =-8; 

RANGELIMIT_UL =40; 


ACTUATOR_N =2; 

ACTUATOR_NUM =[0 0 1]; 

ACTUATOR_DEN =[.0254 .223 1] ; 


VECTORUNITY_N =1; 

VECTORUNITY_U =1 ; 

XCPC_N = 1; 

XCPC_K =10; 


PARENTG 

GOTO FlightControls 
GOTO Collective 

THETLIMIT_N =1; 

THETLIMIT_LL =&WORLD_MODEL_CONTROL_DATA_THETMRLL ; 

THETLIMIT_UL =&WORLD_MODEL_CONTROL_DATA_THETMRUL ; 

TH0CHK_N =1; 

TH0CHK_U =&WORLD__MODEL_CONTROL_DATA_THOCHK ; 

SUMJ4_N =1; 

SUMJ4_SA1 =1; 

SUMJ4_SA2 =1; 


THETA0OUT_N =1; 


BIAS_N =1; 

BIAS_K =&WORLD_MODEL_CONTROL_DATA_THETMRBIAS ; 


SUMJ3_N =1; 
SUMJ3_SA1 =1; 
SUMJ3_SA2 cl; 


KXCTH0_N =1; 

KXCTH0_K =&WORLD_MODEL_CONTROL_DATA_KXCTHO; 
XCILS _N =1; 

XCILS _U = &WORLD_MODEL_CONTROL_DATA_XCS ; 


SUMJ2_N =1; 
SUMJ2_SA1 =1; 
SUMJ2_SA2 =1; 


SUMJ1_N =1; 

SUMJ1_SA1 =1; 
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_SA2 


XC_N =1; 

XC_U =&WORLD_MODEL_CONTROL_DATA_XC; 

XCTRM_N =1; 

XCTRM_U = &WORLD_MODEL_CONTROL_DATA_XCTRM; 


DEG2RAD_N = 1; 

DEG2RAD_K =&W0RLD_DATA_D2R; 


PRISERVO_N =2; 

PRISERVO_NUM = &WORLD_MODEL_CONTROL_DATA_PRIS ERVONUM; 

PR IS ERVO_DEN = &WORLD_MODEL_CONTROL_DATA_PRI S ERVODEN ; 


PARENTG 

GOTO Longitudinal 

SUMJ4_N =1; 
SUMJ 4_SA1 =1; 
SUMJ4_SA2 =1; 

SUMJ 8_N =1; 
SUMJ8_SA1 =1; 
SUMJ8_SA2 =1; 


B1SCHK_N =1; 

B1 SCHK_U = &WORLD_MODEL_CONTROL_DATA_B 1SCHK ; 


SUMJ6_N =1 
SUMJ6_SA1 =1 
SUMJ6_SA2 =1 


XB_N =1; 

XB_U =&WORLD_MODEL_CONTROL_DATA_XB; 


SUMJ1_N =1; 

SUMJ1_SA1 =1; 

SUMJ1_SA2 =1; 

B1SLIMIT_N =1; 

B1SLIMIT__LL = &WORLD_MODEL_CONTROL_DATA_B 1SLL; 
B1SLIMIT_UL =&WORLD_MODEL_CONTROL_DATA_BlSUL; 

SUMJ3_N =1; 

SUMJ3_SA1 =1; 

SUMJ3_SA2 =1| 


KXBB1S_N =1; 

KXBB1S_K =&WORLD_MODEL_CONTROL_DATA_KXBBlS ; 
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PASUM2_N = 1; 

PASUM2_SA1 =1; 
PASUM2_SA2 = -l ; 



=&W0RLD_M0DEL_C0NTR0L_DATA_B1SBIAS; 


=&WORLD_MODEL_CONTROL_DATA_DSAS; 








=&WORLD_MODEL_CONTROL_DATA_PAWOUTNUM; 
=&WORLD_MODEL_CONTROL_DATA_PAWOUTDEN; 


RLD_MODEL_CONTROL_DATA_ASAS; 


RLD_MODEL_CONTROL_DATA_DSAS; 


RLD_MODEL_CONTROL_DATA_PDLLIMIT ? 
RLD_MODEL_CONTROL_DATA_PDULIMIT; 


RLD_MODEL_CONTROL_DATA_PDLAGNUM; 
RLD_MODEL_CONTROL_DATA_PDLAGDEN; 


RLD_MODEL_CONTROL_DATA_PDGAINl; 

RLD_M0DEL_C0NTR0L_DATA_PDGAIN2; 

=&WORLD_MODEL_CONTROL_DATA_PDWOUTNUM; 
= &WORLD_MODEL_CONTROL_DATA_PDWOUTDEN; 
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DEG2 IN_K 


=&W0RLD_M0DEL_C0NTR0L_DATA_KXBB1 SINV; 


PER2DEG_N =1 ; 

PER2DEG_K =&WORLD_MODEL_CONTROL_DATA_KPCDGP; 
PSASLIMIT_N =1; 

PSASLIMIT_LL =&WORLD__MODEL_CONTROL_DATA_PSASLLIMIT ; 

PSASLIMIT_UL =&WORLD_MODEL_CONTROL_DATA_PSASULIMIT ; 


PSASACTUATOR_N =2; 

PSASACTUATOR_NUM =&WORLD_MODEL_CONTROL_DATA_PSASACTNUM ; 
PSASACTUATOR_DEN =&WORLD_MODEL_CONTROL_DATA_PSASACTDEN ; 


PARENTG 

SUMJ2_N =1; 

SUMJ2_SA1 =1; 

SUMJ2_SA2 =1; 

GOTO PBA 

ZERO_N =1; 

ZERO_U =0; 

PBASWITCH_N =1; 

PBASWITCH_SW =&WORLD_MODEL_CONTROL_DATA_PBASW ; 

PBABIAS_N =1; 

PBABIAS_K =&WORLD_MODEL_CONTROL_DATA_PBABIAS ; 


PBATHETLIMIT_N =1 ; 

PBATHETLIMIT_LL = &WORLD_MODEL_CONTROL_DATA_PB ATH ETL LIMIT ; 
PBATHETLIMIT_UL =tWORLD_MODEL_CONTROL_DATA_PBATHETULIMIT ; 


PBAQLIMIT_N =1; 

PBAQLIMIT_LL =&WORLD_MODEL_CONTROL_DATA_PBAQLLIMIT ; 

PBAQLIMIT_UL =&WORLD_MODEL_CONTROL_DATA_PBAQULIMIT ; 

PER2DEGS_N =1; 

PER2DEGS_K =&WORLD_MODEL_CONTROL_DATA_KPCDGP ; 


PBALIMXT2_N =1; 

PBALIMIT2_LL 

PBALIMIT2_UL 

PBALIMIT1_N =1; 

PBALIMIT1_LL 

PBALXMIT1_UL 

INTEGRATOR_N 


=&WORLD_MODEL_CONTROL_DATA_PBALLIMIT2; 
=&WORLD_MODEL_CONTROL_DATA_PBAULIMIT2; 


= &WORLD_MODEL_CONTROL_DATA_PBALLIMITl ; 
=&WORLD_MODEL_CONTROL_DATA_PBAULIMITl ; 
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PRISERVO_N =2; 

PRIS ERVO_NUM =&WORLD_MODEL_CONTROL_DATA_PRI S ERVONUM; 

PRISERVO_DEN =&WORLD_MODEL_CONTROL_DATA_PRISERVODEN ; 


PARENTG 
GOTO Lateral 

A1SCHK_N =1; 

A1SCHK_U =&WORLD_MODEL_CONTROL_DATA_AlSCHK ; 

SUMJ5_N =1; 

SUMJ5_SA1 =1; 

SUMJ5_SA2 =1; 

SUMJ1_N =1; 

SUMJ1_SA1 =1; 

SUMJ1_SA2 = 1; 

XATRM_N =1; 

XATRM_U =&WORLD_MODEL_CONTROL_DATA_XATRM ; 

XA_N =1; 

XA_U = &WORLD_MODEL_CONTROL_DATA_XA ; 


A1SLIMIT_N =1; 
A1SLIMIT_LL =0; 
A1SLIMIT_UL =10.0; 


SUMJ2_N =1; 
SUMJ2_SA1 =1; 
SUMJ2_SA2 =1; 


KXAA1S_N =1; 

KXAA1S_K =&WORLD_MODEL_CONTROL_DATA_KXAAlS ; 


SUM J 4_N =1; 
SUMJ4_SA1 =1; 
SUMJ4_SA2 =1; 


BIAS_N =1; 

BIAS_K =&WORLD_MODEL_CONTROL_DATA_AlSBIAS ; 

AlSOUT_N =1; 


DEG2RAD_N =1; 

DEG2RAD_K =&WORLD_DATA_D2R ; 


GOTO ROLLANALOGSAS 

RAGROUND_N =1 ; 
RAG ROUND_K =0; 
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RASUM3_N =1; 

RASUM3_SA1 =1; 

RASUM3_SA2 =-l; 

ASAS_N =1 ; 

ASAS_SW =&WORLD_MODEL_CONTROL_DATA_ASAS; 

DSAS_N =1 s 

DSAS_SW =&WORLD_MODEL_CONTROL_DATA_DSAS; 

RALIMIT3_N =1; 

RALIMIT3_LL =&W0RLD_M0DEL_C0NTR0L_DATA_RALLIMIT3 ; 
RALIMIT3_UL =&W0RLD_M0DEL_C0NTR0L_DATA_RAULIMIT3 

RALAG_N =1; 

RALAG_NUM =&WORLD_MODEL_CONTROL_DATA_RALAGNUM ; 
RALAG_DEN =&WORLD_MODEL_CONTROL_DATA_RALAGDEN ; 

RASUM1_N =1; 

RAS0M1_SA1 = 1; 

RASUM1_SA2 =1; 

RAGAIN3_N =1; 

RAGAIN3_K =&W0RLD_M0DEL_C0NTR0L_DATA_RAGAXN3 ; 

RASUM2_N =1; 

RASUM2_SA1 =1; 

RASUM2_SA2 =1; 

RAGAIN2_N =1; 

RAGAIN2_K =&W0RLD_M0DEL_C0NTR0L_DATA_RAGAIN2 ; 
RALIMIT2_N =1; 

RALIMIT2_LL =&WORLD_MODEL_CONTROL_DATA_RALLIMIT2 ; 
RALIMIT2_UL =&W0RLD_M0DEL_C0NTR0L_DATA_RAULIMIT2; 

RAGAIN1_N =1; 

RAGAIN1_K =&WORLD_MODEL_CONTROL_DATA_RAGAIN1 ; 
RALAG2_N =1; 

RALAG2_NUM =&WORLD_MODEL_CONTROL_DATA_RALAG2NUM ; 
RALAG2_DEN =&WORLD_MODEL_CONTROL_DATA_RALAG2DEN ; 

RALIMIT1_N =1; 

RALIMIT1_LL =&WORLD_MODEL_CONTROL_DATA_RALLIMITl; 
RALIMIT1_UL =&WORLD_MODEL_CONTROL_DATA_RAULIMIT1 ; 
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PARENTG 

GOTO ROLLDIGITALSAS 

RDGROUND_N =1; 
RDGROUND_K =0 ; 

RDGAIN4_N =1; 
RDGAIN4_K =2; 

RDSUM3_N =1; 
RDSUM3_SA1 =1; 
RDSUM3_SA2 =-l; 


DSAS_N =1; 

DSAS_SW =&WORLD_MODEL_CONTROL_DATA_DSAS ; 


ASAS_N =1; 

ASAS_SW =&WORLD_MODEL_CONTROL_DATA_ASAS ; 


RDLIMIT2_N =1; 

RDLIMIT2_LL =&WORLD_MODEL_CONTROL_DATA_RDLLIMIT2 
RDL IMIT2_UL = &WORLD_MODEL_CONTROL_DATA_RDUL IMIT2 

RDSUM2_N el; 

RDSUM2_SA1 =1; 

RDSUM2_SA2 =1; 


RDGAIN3_N =1 ; 

RDGAIN3_K =&WORLD_MODEL_CONTROL_DATA_RDGAIN3 ; 


RDSUM1_N =1; 

RDSUM1_SA1 =1; 

EDSUM1_SA2 =1; 

RALAG_N =1; 

RALAG_NUM = [ 0 1 ] ; 

RALAG_DEN =[1 1] ; 

RDLIMIT1_N =1; 

RDLIMIT1_LL =&WORLD_MODEL_CONTROL_DATA_RDLLIMIT1 
RDLIMIT1_UL =&WORLD_MODEL_CONTROL_DATA_RDULIMXTl 


RDGAIN2_N =1; 

RDGAIN2_K = &WORLD_MODEL_CONTROL_DATA_RDGAIN2 ; 


RDGAIN1_N =1 ; 

RDGAIN1_K =&WORLD_MODEL_CONTROL_DATA_RDGAIN1 ; 
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LD_M0DEL_C0NTR0L_DATA_KXAA1 SINV ; 

LD_MODEL_CONTROL_DATA_KPCDGR ; 

=&WORLD_MODEL_CONTROL_DATA_RSASLLIMIT; 

=&WORLD_MODEL_CONTROL_DATA_RSASULIMIT; 

=2; 

=&WORLD_MODEL_CONTROL_DATA_RSASACTNUM; 

=&WORLD_MODEL_CONTROL_DATA_RSASACTDEN; 


=&WORLD_MODEL_CONTROL_DATA_PRISERVONUM; 
=&WORLD_MODEL_CONTROL_DATA_PR IS ERVODEN ; 


jD_MODEL_CONTROL_DATA_KXPTHRINV; 

J D_MODEL_CONTROL_DATA_KPCDGY; 

= 1 ; 

=-10.0; 

=10.0; 
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YSASACTUATOR_N =2; 

YSASACTUATOR_NUM =&WORLD_MODEL_CONTROL_DATA_YSASACTNUM 
Y S AS ACTUATOR_DEN =&WORLD_MODEL_CONTROL_DATA_YSASACTDEN 


PARENTG 

GOTO YAWDIGITALSAS 

YDSUM4_N = 1; 
YDSUM4_SA1 =1; 
YDSUM4_SA2 =1; 


YDWASHOUT_N =1; 

YDWASHOUT_NUM =&WORLD_MODEL_CONTROL_DATA_YDWOUTNUM; 

YDWASHOUT_DEN =&WORLD_MODEL_CONTROL_DATA_YDWOUTDEN ; 


YDGAIN3_N =1; 

YDGAIN3_K =&WORLD_MODEL_CONTROL_DATA_YDGAIN3; 


YDLAG2_N =1; 

YDLAG2_NUM =&WORLD_MODEL_CONTROL_DATA_YDLAG2NUM ; 
YDLAG2_DEN =&WORLD_MODEL_CONTROL_DATA_YDLAG2 DEN; 


YDGAIN4_N =1; 

YDGAIN4_K =&WORLD_MODEL_CONTROL_DATA_YDGAIN4; 


YDGAIN2_N =1; 

YDGAIN2_K =&WORLD_MODEL_CONTROL_DATA_YDGAIN2 ; 


ZERO_N =1; 

ZERO_U =0; 


VX5 0KTSWITCH1_N = 1; 

VX50KTSWITCH1_SW =&WORLD_MODEL_CONTROL_DATA_V50KTS ; 


YDLAG1_N si; 

YDLAG1_NUM =&WORLD_MODEL_CONTROL_DATA_YDLAGlNUM ; 
YDLAG1_DEN =&WORLD_MODEL_CONTROL_DATA_YDLAG1DEN ; 


YDSUM1_N =1; 
YDSUM1_SA1 = 1; 
YDSUM1_SA2 =1; 


YDGAIN1_N =1; 

YDGAIN1_K =&WORLD_MODEL_CONTROL_DATA_YDGAIN1 ; 


ASAS_N =1; 

ASAS_SW -&WORLD_MODEL_CONTROL_DATA_ASAS ; 

DSAS_N =1; 

DSAS_SW =&WORLD_MODEL_CONTROL_DATA_DSAS ; 
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YDSUM3_N = 1; 
YDSUM3_SA1 =1; 
YDSUM3_SA2 =-l; 

YDGAIN5_N =1; 
YDGAIN5_K =2 ; 

YDGROUND_N =1 ; 
YDGROUND_K =0; 

YDSUM2_N =1; 
YDSUM2_SA1 =1; 
YDSUM2_SA2 =1; 

YDUNITY_N =1; 
YDUNITY_U =1; 

YAWSWITCH_N =1; 

YDSUM6_N =1; 

YDSUM6_SA1 =-l; 
YDSUM6_SA2 =1; 

YDGAIN_N =1; 
YDGAIN_K =49.0; 

VXKNOTS_N =1; 
VXKNOTS_K =1; 


YDLIMIT_LL =&WORLD_MODEL_CONTROL_DATA_YDLLIMIT; 
YDLIMIT_OL =&WORLD_MODEL_CONTROL_DATA_YDULIMIT; 


PARENTG 

THRCHK_N 

THRCHK_U 

SUMJ5_N 

SUMJ5_SA1 

SUMJ5_SA2 

SUMJ1_N 

SUMJ1_SA1 

SUMJ1_SA2 


=&WORLD_MODEL_CONTROL_DATA_THRCHK; 


=&WORLD_MODEL_CONTROL_DATA_XPTRM; 


I 


I 




XP_U =&WORLD_MODEL_CONTROL_DATA_XP; 


THETTRLIMIT_N 
THETTRL IMIT_LL 
THETTRLIMIT_UL 


=&WORLD_MODEL_CONTROL_DATA_THETTRLL ; 
=&WORLD_MODEL_CONTROL_DATA_THETTRUL; 


SUMJ2_N =1; 
SUMJ2_SA1 =1; 
SUMJ2_SA2 =1; 


KXPTHR_N =1; 

KXPTHR_K = &WORLD_MODEL_CONTROL_DATA_KXPTHR; 


SUMJ 4_N =1; 
SUMJ 4_SA1 =1; 
SUMJ4_SA2 =1; 


BIAS_N =1; 

B IAS_K =&WORLD_MODEL_CONTROL_DATA_THETTRBI AS ; 


THETTROUT_N =1; 


DEG2RAD_N =1; 

DEG2RAD_K =&W0RLD_DATA_D2 R; 


GOTO YAWANALOGSAS 

YAWASHOUT_N =1; 

YAWASHOUT_NUM = &WORLD_MODEL_CONTROL_DATA_YAWOUTNUM 

YAWASHOUT_DEN =&WORLD_MODEL_CONTROL_DATA_YAWOUTDEN 


ZERO_N =1; 

ZERO_U =0; 


YAGAIN4_N =1; 

YAGAIN4_K =&WORLD_MODEL_CONTROL_DATA_YAGAIN4; 


YASUM2_N =1 
YASUM2_SA1 =1 
YASUM2_SA2 =1 


VX50KTSWITCH2_N =1; 

VX5 0KTSWITCH2_SW =&WORLD_MODEL_CONTROL_DATA_V5 OKTS; 
YAGAIN3_N =1 ; 

YAGAIN3_K = &WORLD_MODEL_CONTROL„DATA_YAGAIN3 ; 


VX5 0KTSWITCH1_N =1; 

VXS0KTSWITCH1_SW =&WORLD_MODEL_CONTROL_DATA_V5 OKTS ; 
YALIMIT3_N =1; 

YALIMIT3_LL =&WORLD_MODEL_CONTROL_DATA_YALLIMIT3 ; 
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YAGAIN5_N =1; 
YAGAIN5_K =2 ; 

YAGROUND_N =1; 
YAGROUND_K =0; 

YASUM1_N = 1; 
YASUM1_SA1 =1; 
YASUM1_SA2 =1; 


YAGAIN2_K =&WORLD_MODEL_CONTROL_DATA_YAGAIN2 ; 
YAGAIN1_N =1; 

YAGAIN1_K =&WORLD_MODEL_CONTROL_DATA_YAGAIN1; 
YALAG_N =1; 

YALAG_NUM =&WORLD_MODEL_CONTROL_DATA_YALAGNUM; 
YALAG_DEN = &WORLD_MODEL_CONTROL_DATA_YALAGDEN ; 

YALIMIT1_N =1; 

YALIMIT1_LL =&WORLD_MODEL_CONTROL_DATA_YALLIMITl; 
YALIMIT1_UL =&WORLD_MODEL_CONTROL_DATA_YAULIMITl ; 

YALIMIT2_N =1; 

YALIMIT2_LL =&WORLD_MODEL_CONTROL_DATA_YALLIMIT2 ; 
YALIMIT2_UL =&WORLD_MODEL_CONTROL_DATA_YAULIMIT2 ; 


PARENTG 





PR IS ERVO_NUM =&WORLD_MODEL_CONTROL_DATA_PRISERVONUM ; 

PRI SERVO__DEN =&WORLD_MODEL_CONTROL_DATA_PRISERVODEN ; 

PARENTG 

PARENTG 

PARENTG 

PARENTG 
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world_data_xa 
world_da ta_xb 
world_data_xc= 
world_daca_xp=&xp, 

Xatrm=4.95; 
Xbtrro=l.196; 
Xctrm=5.00; 

Xptrm=l.925; 
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// yaw digital inputs 
ydwoutnum=[2 0]; 
ydwoutden=[2 1]; 
ydlaglnum=[0 1]; 
ydlaglden=[l.045 1] ; 
ydlag2num=[0 1]; 
ydlag2den=[1.045 1]; 
ydgainl=-.611; 
ydgain2=-l.013; 
ydgain3=-l.23; 
ydgain4=.152; 
ydllimit=-5.0; 
ydulimit=5.0; 

// yaw sas input to mechanical controls 
ysasactnum=[0 0 1]; 
ysasactden=[.00032 .025 1]; 
ysasllimit=-10.0; 
ysasulimit=10.0; 

kpcdgy=-.298; 
kxpthrinv=l/kxpthr; 

// Inputs for Pitch Bias Actuator 
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APPENDIX E. OPEN LOOP EIGENANALYSIS 


x= [ u ; w; q; 0; V; p; <J); r ] 
u=[6.;6,;6 c ;8 p ] 

Equivalent Airspeed = 0 kts, Gross Weight = 19462 lbs, SSL Conditions 
A= 

[0.2032 0.0140 1.1439 -32.1410 -0.2340 -1.7236 -0.0001 -0.0341; 
0.2417 -0.2776 -0.0737 -1.9453 -0.2371 -0.4255 1.4656 2.7489; 










0.0031i 0.0020 + 0.003H 0.0370 


Columns 1 through 4 


-2.9590 

0 

0 


0 0 
0 0 


Columns 5 through 8 


0 0 

0.4479 + 0.8051i 0 

0 0.4479 - 0.80511 

0 0 


0 0 
0 0 


0 

0 


0.7587i 0 0 

0.0073 - 0.7587i 0 

0 -0.1569 

0 0 


0 


-0.2659 
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Columns 1 through 4 

-2.8769 0 

0 -1.2267 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

Columns 5 through 8 


0 0 

0 0 

-0.1584 + 0.94401 0 

0 -0.1584 

0 0 

0 0 

0 0 

0 0 


0.9440i 


0 

0 

0 

0 

0.4249 

0 




0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

-0.2667 + 0.3470i 0 0 

0 -0.2667 - 0.3470i 0 


0 0 -0.1539 
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eigval = 

Columns 1 through 

-3.5311 

0 


0 

0 

0 

0 

0 
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